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ABSTRACT 

We present optical and near-infrared adaptive optics (AO) imaging and spectroscopy 
of 13 ultracool (>M6) companions to late-type stars (K7-M4.5), most of which have 
recently been identified as candidate members of nearby young moving groups (YMGs; 

8-120 Myr) in the literature. Three of these are new companions identified in our AO 
imaging survey and two others are confirmed to be comoving with their host stars 
for the first time. The inferred masses of the companions (~ 10-100 Mjup) are highly 
sensitive to the ages of the primary stars so we critically examine the kinematic and 
spectroscopic properties of each system to distinguish bona fide YMG members from old 
field interlopers. 2MASS J02155892-0929121 C is a new M7 substellar companion (40- 
60 Mjup) with clear spectroscopic signs of low gravity and hence youth. The primary, 
possibly a member of the ~40 Myr Tuc-Hor moving group, is visually resolved into three 
components, making it a young low-mass quadruple system in a compact (<100 AU) 
configuration. In addition, Li I A6708 absorption in the intermediate-gravity M7.5 
companion 2MASS J15594729-I-4403595 B provides unambiguous evidence that it is 
young (<200 Myr) and resides below the hydrogen burning limit. Three new close- 
separation (<1^0 companions (2MASS J06475229-2523304 B, PYC J11519-I-0731 B, 
and GJ 4378 Ab) orbit stars previously reported as candidate YMG members, but 
instead are likely old (>1 Gyr) tidally-locked spectroscopic binaries without convincing 
kinematic associations with any known moving group. The high rate of false positives 
in the form of old active stars with YMG-like kinematics underscores the importance of 
radial velocity and parallax measurements to validate candidate young stars identified 
via proper motion and activity selection alone. Finally, we spectroscopically confirm 
the cool temperature and substellar nature of HD 23514 B, a recently discovered M8 
benchmark brown dwarf orbiting the dustiest-known member of the Pleiades. 

Subject headings: binaries; spectroscopic — brown dwarfs — stars: individual (2MASS J02155892- 
0929121, 2MASS 415594729-^4403595, HD 23514) 
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Introduction 


With masses between ~13-75 Mj^p, brown dwarfs directly link the lowest-mass stars and 
young giant planets in mass, radius, and temperature. Nearly two thousand isolated brown dwarfs 
have been identified over the past 20 years, but their individual ages, metallicities, and masses are 
largely unknown. Age benchmarks — brown dwarfs orbiting well-characterized stars or members of 
coeval clusters — are rare but valuable tools to test low-temperature atmospheric and evolutionary 
models for mutual consistency and empirically calibrate isolated field objects b ecause their ages , 
distances, and compositions can be le veraged from the host star or cluster (e.g.. lNaud et al.ll2ni4l : 
Crepp et al. 2014 : Deacon et al. 2014 ). 


Over the past decade, dozens of nearby brown dwarfs with optical and near-infrared (NIR) 


spectroscopic signs of l ow surface gravity (youth) have been identified (e.g., iRebolo et al. 


1995 


Kirkpatrick et al 


2008 


Allers &: Liull2013l i. Many of these objects have proper motio ns consistent 


with nearby (<80 pc) young moving groups (YMGs) with ages of ~10-120 Myr (e.g. . IScholz et al 


2OO5I : Ihooper et ahl 120071 : IPelorme et al.ll2012l : iFahertv et al.ll2012l : lOagne et al.ll2014l ) but because 


they are intrinsically faint, most of these isolated brown dwarfs lack trigonometric distances and 
radial velocities (RVs) needed to unambiguously define their space motion. As a result, only 


YMGs like TWA 

Gizis et al. 

to 

0 

0 

Weinbereer et al. 

2013), B Pic ( 

Liu et al. 

2013b), and AB 

Dor ( 

Fahertv et al. 

2 OI 2 I: 

Liu et al. 2 

013a 

). These low-gravity benchmarks are especially useful for 


empirical spectral classification sequences to map the effects of age on gravity-sensitive features like 
alkali absorption lines, i7-band morphology, and H 2 O steam bands at ~1.4 and ~1.9 /rm spanning 
a large range of temperatures (500 K-2500 K) and corresponding masses (^^5-75 Mjup). These 
empirically-based spectral sequences at hxed ages also enable detailed analyses of secondary effects 



Another approach to find these elusive young substellar benchmarks is to search for faint 
companions to known nearby young stars. One advantage of this method is that precise paral¬ 
laxes and RVs are more readily obtained for the primary stars, making kinematic association with 
YMGs relatively straightforward compared to faint isolated brown dwarfs in the held. Indeed, 


space-based imaeine (Lowrance et al. 

. 999 I: Webb et al. 

I 999 I: Lowrance et al. 

200011 and adaptive 

optics (AO) imaging from the ground (< 

Fhauvin et al. 

2 OO. 5 I: E 

Idler et al. 

20ld 

Mugrauer et al. 201C: 

Wahhai et al. 2011: Bowler et al. 2012c 

j: Delorme et al. 

2013 

: Bowler et al. 2 

OI 3 IL However, old in- 


terloping stars with UVW space motions similar to YMGs can complicate efforts to identify young 
substellar companions in moving groups. Distinguishing bona hde members from kinematic in¬ 
terlopers can be particularly challenging for M-type members of YMGs since both single mid- to 
late-M dwarfs and tidally locked M-|-M binaries can have short rotation periods (less than a few 
days) and strong magnetic dynamos even at old (^1 Gyr) ages. 
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Here we present a suite of spectroscopic and adaptive optics imaging observations of ultracool 
companions to M dwarfs. The majority of these companions have been identified as part of the 
Planets Around Low-Mass Stars (PALMS) high-contrast imaging search for giant planets and brown 
dwarfs, three of which are new here (2MASS J06475229-2523304 B, PYC J11519-I-0731 B, and 
GJ 4378 Ab). Seve ral additional systems were identified by the Astralux Lucky imaging survey 
(jjanson et al.l 120121 ) . which obtained shallow diffraction-limited imaging of nearby, mostly active 
M dwarfs in the optical. The sky positions, proper motions, spectral types, distances, and ages 
of our targets are listed in Table [TJ The physical properties of the companions can be found in 
Table [21 Most of the primary stars have been previously reported in the literature as candidate or 
confirmed members of YMGs, making the companions excellent candidates for young benchmark 
brown dwarfs. 


We critically examine these moving group membership claims on a case-by-case basis using 
new RVs of the primary stars together with recently published parallaxes, when available. We also 
use our new spectroscopy of the companions to independently constrain the ages of the systems 
via low-gravity spectroscopic features, which can be more prominent and last longer (<200 Myr) 
at late-M and L spectral types than at earlier types. Finally, we search for signs of orbital motion 
with multi-epoch high-resolution AO imaging to identify systems that will yield dynamical masses 
in the near future. 


2. Observations 


2.1. Adaptive Optics Imaging 

2.1.1. Keck/NIRC2 


Our near-infrared imaging with NIRG2 at the 10-meter Keck II teles cope was carried gut be - 
tween 2012 and 2014 using Natural Guide Star Adaptive Optics (NGSAO; Wizinowich et al. 2000). 
We utilized the na rrow NIRG2 camera mode, resulting in a plate scale of 9.952 ± 0.002 mas pix“^ 
( Yelda et al. 201C 1 and a field of view of lC/.'2 x lCf! 2. The NIRG2 J, H, K, Ko, and L' filters are 
on th e Maunakea Observatory (MKO) filter system ([Simons &: Tokunagall2002l : iTokunaea &: Vacca 
2005). The NIRC2 Y-band filter was installed in 2011 and has a central wavelength of 1.0180 Min 
wit h a bandwidth o f 0.0996 /rm, differing slightly from other Y-band filters ([Hillenbrand et al.ll2002l : 
see iLiu et al.ll2012l for additional details). In addition, we also include a single epoch of A'5-band 


NGSAO imaging of IRXS J034231.8+121622 from 2007 in our analysis of the system ('Section 13.41) . 

Raw images were cleaned of cosmic rays and bad p ixels, dark-subtrac ted, flat-fielded, and 
north-aligned using the orientation of +0?252 ± 0?009 from I Yelda et al.l (|2010l L For our L' data, the 
science frames were combined after masking the target to create a flat field. For close companions 
(<2"), astrometry and relative pho tometry are mea sured by fitting three bivariate Gaussians to each 
binary component as described in iLiu et al.l (|2008l L Aperture photometry is used for companions 





























widely-separated from their host stars. The standard deviations from multiple measurements are 
adopted for our separation and position angle (P-A.) uncertainties. 


For PYC J11519+0731 we also ob tained deep imaging with the 6 00-mas diameter coronagraph 
in angular differential imaging mode (jLiul 120041 : iMarois et al.l 120061 ) in H band on 2012 May 22 
UT. We acquired 20 frames each with an integration time of 30 sec spanning a total sky rota¬ 
tion of 11°. Following image registration based on the star’s position behind the partly opaque 
mask and corrections for optical distortions (B. Cameron, 2007, priv. comm.), w e performed PSF 


subtr action with the Locally-Optim i zed C ombination of Images algorithm (LOCI; iLafreniere et al 


20071 1 as described in iBowler et al.l (|2015lL However, because of the small sky rotation, we also 


used a reference PSF library containing over 2000 NIRC2 coronagraphic images from our PALMS 
program. At each geometric subsection the LOCI coefficients were computed from 100 library im¬ 
ages, which were chosen by minimizing the rms after scaling and subtracting every PSF reference 
image within that particular section. For our final reduction we masked the region surrounding 
PYC J11519-I-0731 B to avoid s elf-subtractio n when computing image coefficients. Coronagraph 
throughput measurements from iBowler et al.l (j2015|) were used to produce a 7 a contrast curve. 
Results are described in Section 13.111 


2.1.2. Subaru/IRCS 


We imaged 2MASS J02155892-0929121 2MASS J06475229-2523304, and GJ 4378 A with the 


Infrared Camera and Spectrogr aph (IRCS; 


Tokunaga et al. 


1998 


Kobavashi et al.ll200d l coupled 


with the A0188 NGSAO system ( Havano et al.ll20iri ) at the 8.2-meter Subaru Telescope on UT 2012 
October 2012 12 and 13. The weather was clear throughout both nights but the seeing was highly 
variable (1-2''). The A-B-C components of 2MASS J02155892-0929121 were easily identified, but 
the (then unknown) fourth member, Ab, was not resolved in these data. We obtained MKO iL-band 
images of 2MASS J02155892-0929121 ABC and MKO J, H, and K images of 2MASS J06475229- 
2523304 AB and GJ 4378 Aab. 

The raw frames were cleaned of bad pixels, flat-fielded, and dark-subtracted. Astrometry and 
relative photometry were measured f or ea ch image individually by fitting three bivariate Gaus- 


sians to each component (iLiu et al.l 120081 ). We adopt a plate scal e and North orientation of 


20.41 ± 0.05 mas pix“^ and 4 -89?03 ± 0?09 measu red in IBowler et al.l ( 20131 ) using the 2'.'4 binary 
IRXS J235133.3-1-312720 AB ( Bowler et al. 2012a ). which we observed with IRCS on UT 2012 Oc¬ 
tober 12. 


2.1.3. Gemini/NICI 


On UT 2010 Aug ust 29 we targeted IRXS J034231.8-1-121622 with the Near-Infrared Corona- 
graphic Imager (NICI; IChun et al.ll2008l l on the 8.1-meter Gemini South telescope in queue mode 
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as part of the NICI Planet-Finding Campaign (jLiu et alJl2010bl : Program GS-2010B-Q-500, PI: M. 
Liu). NICI is a simultaneous dual-channel imager with a tapered, partly opaque Lyot coronagraph 
and two 1024x1024 (18"xl8" field of view) infrared arrays. A short imaging sequence of eight 
60 s frames was acquired for IRXS J034231.8-1-121622 in H and CH^long (4% throughput cen¬ 
tered at 1.652 /im) bandpasses. Ultimately this target it was not included in the final NICI Planet 
Finding Campaign because the ADI sequence was aborted after a f ew frames. The Cf.'8 companion 
IRXS J034231.8-1-121622 B originally found by Bowler et al. ( 2015 ) is easily detected in individual 
images without PSF subtraction. 


PYC J11519-I-0731 was targeted in a similar fashion with NICI on 2013 April 19 UT as part 
of a separate survey to search for giant planets around low-mass stars (Program GS-2013A-Q- 
54, PL E. Nielsen). A total of 10 short (1.14 s per coadd x 10 coadds) frames were obtained 
simultaneously in H and Ks bands, resulting in a total on-source integration time of 114 s per 
filter. PYC J11519-I-0731 B was easily resolved and remained unsaturated in these short images. 


We reduced the NICI images following the description in IWahhai et al.l ([20111) , which includes 
basic image reduction (bad pixel removal and flat fielding) and distortion correction. Plate scale and 
sky orientation values were used fro m calibration measurem ents by the Gemini staff closest in time 
before the observations were taken (jHavward et al.M2014l ). Astrometry a.nd relative photome try 
incorporating the mask transmission measurement and uncertainty from IWahhaj et al.l (|201lh is 
listed in Table O 


2.1.4- Palomar 60-Inch/Robo-AO 


2MASS J15594729-h4403595, 2MASS 411240434-^3808108, and G 180- 11 were imaged with 


the visible-light Ro bo-AO robotic laser guide star adaptive optics system (IBaranec et al 


2013; 


Baranec et al.l[2014 l mounted on the automated Palomar 60-inch (1.5-meter) telescope ( Genko et al 


2003) on the night of 2014 June 13 UT. Robo-AO utilizes a pulsed ultraviolet laser to generate a 
Rayleigh-scattered guide star for wavefront sensing, enabling efficient observations with AO correc¬ 
tion in the optical down to U~16 mag. The electron-multiplying GGD camera is conti nually read 


out d uring the observations and a shift-and-add pipeline creates a final science im age (ILaw et al 
2014). The nominal detector field of view is 44" with a plate scale Cf.T431 pixeU^ ( Baranec et al 


20141 ). The pipeline processing produces a final image that is resampled at half the original pixel 
scale (0'.'0216 pixeU^). 


All three targets were imaged in Sloan r', i', and z' filters ([York et al.ll2000l ) with exposure times 
of 60 sec each. The late-M companions to 2MASS J15594729-I-4403595 and 2MASS J11240434-I-3808108 
are clearly visible in the i' and z' images. G 180-11 B is detected in all three filters. The North 
orientation and plate scale are derived us i ng th e z' image of 2MASS J15594729-I-4403595 AB and 
astrometry of the pair from iBowler et al.l (120151 ). Based on its photometric distance of 27 pc, this 
system has a projected separation of ~150 AU. No orbital motion is expected as the observations 











































used for astrometric calibration were taken in 2012. We measure a plate scale of C/.'0219 ± C/.^0002 
pixel“^ and sky orientation on the detector of 335?4 ± 0?4. These are within 2 u of the values de¬ 
rived using a distortion solution based on ff5T-calibrated images of M15 (A. Tokovinin, 2015, priv. 
communication). The FWHM of the primary stars in these observations are ~0l'28, ~C/.'20, and 
~C/.'16 in r', i', and z', respectively, in 1"3-1"8 seeing conditions. Flux ratios for the components 
are measured with aperture photometry and listed in Table [3l 


2.2. Near-Infrared Spectroscopy 

2.2.1. Keck/OSIRIS-NGSA O 


Near-infrared spectroscopy of low-mass companions allows for spectral classification and pro¬ 
vides an independent diagnostic of youth in the form of gravity-sensitive atomic and molecular 
features. We obtained moderate-resolution (i2=(5A/A=3800) near-infrared spectra of the close-in 
companions (<3") PYC J11519-h0731 B, G 180-11 B, HD 23514 B, GJ 4378 Ab, 2MASS J06475229- 
2523304 B, and 2MASS J08540240 -3051366 B with the OH-Suppressing Infrared Imaging Spectro¬ 
graph (OSIBIS: iLarkin et al.ll200(Tl ) at the 10-meter Keck I telescope between 2012 June and 2014 
Dec (Table 1^. Our 2013 and 2 014 observations utilized the newly installed, more efficient grating 
described in iMieda et al.l ( 2014 ). We used the 20 mas pixel~^ plate scale (0'.'32xl"28 field of view) 
for PYG J11519-F0731 B, 2MASS J06475229-2523304 B, and 2MASS J08540240-3051366 B. The 
50 mas pixel”^ scale {(f!8x3'.'2 field of view) was used for the other companions (see Table 0]). All 
observations were taken in ABBA nodded pairs for sky subtraction. Immediately before or after 
our science observations we targeted single AOV standards a similar airmass and sky position for 
telluric absorption measurements. 

Basic 2D image reduction, pairwise sky subtraction, and rectification into 3D cubes was car¬ 
ried out with version 3.2 of the OSIRIS Data Reduction Pipeline maintained by Keck Observa¬ 
tory. Rectification matrices derived closest in time prior to the observations were used for the 
pipeline processing. The science and standard spectra were extracted using aperture photome¬ 
try, normalized to a common level, then median combined. Telluric correction was carried out 


with xtellcor_general, p art of the Spextool reduction package for IRTF/SpeX (jVacca et al 


2003; 


Cushing et al.l 120041 '). Individual bandpasses were flux-calibrated with photometry from 


Rodriguez et al.l (120121 1 for HD 23514 B and from Table [3] for the remaining targets. The me¬ 
dian S/N per pixel of our (unsmoothed) OSIRIS spectra range from 25-120. 


2.2.2. IRTF/SpeX 


We obtained moderate-resolution spectra of the wide-separation {3"-8") companions 2MASS 
J11240434-h3808108 B, G 160-54 B, and 2MASS J02155892-0929121 C with the Infrared Tele- 
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scope Facility’s SpeX spectrograph in SXD mode on the nights of UT 2011 January 22, UT 2012 
September 24, and UT 2014 January 18, respectively. The Cf.'5 slit width was used for all three 
targets, yielding a resolving power of i?~1200. We also observed the unresolved primary star 
PYC J11519+0731 Aab in SXD mode with a C/.'8 slit width (J2~750) on UT 2012 June 1. All 
science data were taken in a nodded ABBA pattern. AOV standards, internal flats, and arc lamps, 
were obtained at a similar sky position to the science target. Individual spect ra were extracted. 


median-combined, a nd telluric-corrected with the Spextool reduction package (IVacca et al 


2003 


Cushing et al.l 120041 ). The median S/N per pixel of our SpeX spectra range from 30-60. 


2.3. Optical Spectroscopy 

2.3.1. Mayall/RC-Spectrograph 


We obtained low-resolution optical spectra of most of the primary stars in our sample to provide 
uniform spectral classification and measure Ha line strengths. We targeted PYC J11519-I-0731 Aab 
with the Ritchey-Chretien Spectrograph (RC-Spec) using the T2KA CCD mounted on the NOAO 
Mayall 4-meter telescope on 2014 May 22 UT. The l"5x98" slit was hxed in the North-South 
direction throughout the night so we observed PYC J11519-I-0731 at a low air mass (sec z = 1.1) 
near transit to minimize slit losses from differential atmospheric refraction (jFilippenkdll982l ). We 
used the 600 1 mm“^ BL420 grating in first order which provided a resolving power of i?~2600 
spanning 6300-9200 A. The raw 240 s image was corrected for bad pixels and cosmic rays, bias- 
subtracted, and flat-fielded. We then extracted the spectrum by summing in the spatial direction 
and corrected for throughput losses using our observations of the spectrophotometric standard 
HZ 44. 


2.3.2. SOAR/Goodman Spectrograph 

On 2014 June 28 UT we obtained low resolution optical spectra of the (unresolved) primary 
stars 2MASS J01225093-2439505, 2MASS J02155892-0929121, GJ 4378 A, and GJ 4379 B with 
the Goodman High Throughput Spectrograph at the 4.1-meter Southern Astrophysical Research 
(SOAR) Telescope located at Gerro Pachon, Ghile. All targets were observed at low airmass (sec z < 
1.2) with clear sky conditions. The C/.'46 long slit was used with a CCD read rate of 400 kHZ and the 
GG455 order blocking filter. We acquired low-resolution (i?~1800) spectra of all four targets with 
the 400 1/mm M2 grating spanning ~5000-9000 A. We also observed 2MASS J02155892-0929121, 
GJ 4378 A, and GJ 4379 B with a moderate-resolution setup (R~5900) using the 1200 1/mm 
M5 grating spanning ~6300-7500 A. Integration times were between 240-300 sec for all of our 
observations and the slit was rotated to the parallactic angle to avoid slit losses from chromatic 
dispersion. The Goodman spectrograph suffers from strong fringing redward of ~7000 A, which can 
change from target to target from mechanical flexure. To calibrate and correct these systematics 
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we obtained arc lamp spectra for wavelength calibration and quartz flats after each target. For the 
400 1/mm grating we used HgAr lamps and for the 1200 1/mm grating we used CuHeAr lamps. 


After basic image reduction entailing bad pixel correction, cosmic ray removal, bias subtraction, 
and flat fielding, we removed night sky lines in the 2D image using the chip regions adjacent to 
the science spectrum in the spatial directi on. Th e spec trum was then horizontally rectified and 
optimally extracted using the method from iHornd (|l986l i. For each target, fringing was corrected 
using the same chip region occupied by the science spectrum on the flat field. By measuring the 
effects of fringing in real time at each telescope position over the same pixels covered by the science 
data, we effectively divided out most of these artifacts in our target spectra. Finally, we correcte d 
for throughput losses with the standard star LTT 6248 (jHamuv et al.lll992l : iHamuv et al.lll994l '). 
which was observed with both grating settings. 


2.3.3. UH 2.2-meter/SNIPS 


We obtained low-resolution optical spectra of the host stars G 160-54, 2MASS J15594 729-I-4403595 , 


and 2 MASS J06475229-2523304 with the SuperNova Integral Field Spectrograph fSNIFS aLantz et al. 


200J) at the University of Hawai’i 2.2-meter telescope on the nights of UT 2012 Jan 29, UT 2012 
Aug 23, and UT 2014 Dec 3. SNIFS is an optical integral field spectrograph providing low-resolution 
(7Z~1300) spectroscopy spanning 3200-11000 A. Several B and A-type stars and white dwarfs were 
targeted throughout the nights for use as spectrophotometric standards (Table 0]). The short and 
long-wavelength spectra were redu ced, rectified into spectral cubes, and ex tracted in real time with 
the SNIFS data reduction pipeline ( Aldering et al. 2006 : Scalzo et al. 20ld ). This entailed flat field¬ 
ing, bias subtraction, wavelength calibration using arc lamps, and spectra l extraction using a PSF 
model. Additional flux calibration and telluric correction is described in iMann et al.l ( 2013 ). To 
summarize, we used spectrophotom etric standards in c onjunction with a model of the atmospheric 
transparency above Maunakea from iButon et al.l (j2012l ) to derive airmass and time dependent cor¬ 
rections to the spectrum. The resultin g flux calibration has been shown to be accurate to 1-2% 
when compared to space-based spectra (jMann et al.ll2013l l. 


2 . 3 . 4 . Keek/HIRES 


We obtained high-resolution spectra for most of the primary stars in our sample to measure 
RVs, search for Li I A6708 absorption, and check for spectroscopic binarity. Fully convective low- 
mass stars and high-mass brown dwarfs burn lithium on characteristic timescales that depend 
on their mass, so the presence or absence of the li thium feature provides an independent age 
constraint for the system (e.g., Chabrier et al. 1996 L We targeted 2MASS J02155892-0929121 
A, 2MASS J02594789- 0913109 A, and P YC J11519-I-0731 Aab with the High Resolution Echelle 
Spectrometer (HIRES; Vogt et al. 1994) on the Keck I 10-m telescope. We used the 0.861" slit 
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with HIRES to give a spectral resolution of A/AA~58000. The detector consists of a mosaic of 
three 2048 x 4096 15-//m pixel CCDs, corresponding to a blue, green and red chip spanning 4900 
- 9300 A. We used the GG475 filter with the red cross-disperser to maximize the throughput near 
the peak of a M dwarf spectral energy distribution. 

Standard reduction was performed on each stellar exposure, i.e. bias-subtraction, flat-fielding, 
and the correction of bad pixels. After optimal extraction, the 1-D spectra were wavelength cali¬ 
brated with a Th/Ar arc taken within an hour of the stellar exposure. The resulting spectra reached 
S/N of 20 - 50 per pixel at 7000 A, depending on brightness. Each night, spectra were also taken 
of the RV standard GJ 908 and an AOV standard star for telluric line correction. 


2.3.5. CFHT/ESPaD OnS 


We acquired high-resolution echelle optical spectra of G160-54 A, 2MASS J15594729-I-4403595 
A, and PYG J1151 9-I-0731 Aab using t he Echelle SpectroPolarimetric Device for the Observation of 
Stars (ESPaDOnS: iDonati et al.ll2006l ) on the Ganada-France-Hawaii 3.6-m telescope. ESPaDOnS 
is fiber fed from the Cassegrain focus to the Coude focus. The fiber image is projected onto a 
Bowen-Walraven sheer at the entrance to the spectrograph. We used the “star-l-sky” mode, which 
records the full spectrum over 40 grating orders onto a 2048 x 4608-pixel CCD detector covering 3700 
to 10400 A with a spectral resolution Ri68000. The data were reduced using Lihre Esprit, a fully 


automated r eduction package provided for the instrument and described in detail by iDonati et al 
( 1997 . 2006 1. RV standards used for these targets were GJ 273, GJ 628 and GJ625. 


2.3.6. Irenee du Pont/Echelle Spectrograph 

Four of our targets were observed with the optical echelle spectrograph mounted on the 2.5- 
m du Pont telescope at Las Campanas Observator}!!]: 2MASS J06475229-2523304 A, 2MASS 
J08540240-3051366 A, GJ 4378 A, and GJ 4379 B. The TEK#5 CCD with a 2k by 2k format 
and 24 micron pixels was employed for observations prior to December 2009. At this time, the 
SITe2K CCD was installed with the same format. In order to reach a spectral resolution of 45000, 
we used the 0.75"by 4" slit, and acquired simultaneous wavelength coverage from 3700 to 8500 A. 
We obtained a ThAr exposure at every pointing to avoid any flexure effects and day-time “milky” 
flats correct for pixel-pixel variations. In addition to a hot standard star for telluric line correction, 
RV standards GJ 273, GJ 433, GJ 908, and GJ 699 were also observed. 

Excluding our HIRES spectrum of PYG J11519-I-0731 A, which is detailed in Section [3.111 
out RV measurements for the HIRES, ESPaDOnS, and du Pont telescopes are carried out in a 


^The instrument manual can be found here: http://www.lco.cl/telescopes-information/irenee-du-pont/instruments/website/echelle-s 
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similar fashion. After correcting for the heliocentric velocity, we cross-correlated each spectral 
order between 7000 and 9000 A o f each stellar spec trum with an RV standard of similar spectral 
type using IRAF’Jl fxcor routine Fitzpatrick (1993). We excluded the Ca II infrared triplet (IRT) 
(the target stars exhibit Ca II emission that is not present in the RV standards) and regions of 
strong telluric absorption in the cross-correlation. The resulting cross-correlation function also 
allowed to quickly identify double- or triple-lined spectroscopic binaries: G160-54 A, GJ 4378 A, 
and PYG J11519-I-0731 A. In nearly all cases we are able to measure RVs to better than 1 km s“^. 


3. 7. ESO-MPG 2.2-meter/FEROS 

As part of the ongoing GASTOFFS survey to identify young, low-mass stars in the solar 
neighborhood (Schlieder et al. 2015, in prep), PYG J11 519-k0731 Aab was observed with the 
Fiberfed Extended Range Optical Spectrograph (FEROS; Kaufer et al. 1999 1 on the ESO-MPG 
2.2-m telescope located at La Silla, Ghile, on UT 22 February 2013. FEROS is a fiber-fed optical 
echelle spectrograph providing a resolution RpsdSOOO from fti3500-9200 A across 39 orders. We 
acquired three exposures each with integration times of 406 s at an airmass of 1.34 which resulted 
in a S/N per pixel of >50 in the red orders of the median-combined spectrum. FEROS contains 
two 2'.'0 optical fibers separated by 2f9 on sky. We observed in ‘object -|- sky’ mode with one 
hber on the science target and the other on sky. Data reduction was performed using the FEROS 
Data Reduction System (DRS) implemented within the ESO-MIDAS software package. The DRS 
performs flat-fielding, background subtraction, bad pixel removal, order extraction (using optimal 
extraction), and wavelength calibration from ThAr lamp lines. The DRS then computes barycentric 
velocity correctionJ§ and rebins the individual orders to the same wavelength scale to merge them 
and produce a continuous spectrum. The calibrations used by the DRS (bias, flats, and a ThAr 
lamp spectrum) are acqui red in the afternoo n prior to the observations. We also observed late-type 
RV templates drawn from IPrato et al. ( 2002 ) during observing runs in December 2011 and October 
2012 for use during cross-correlation analysis (Section 13. 111) . Each night of FEROS observation s 
als o includes observatio ns of at least one RV standard drawn from the lists of iGhubak et al.l ()2012l ) 


or 


Nidever et al.l (|2002l ) for internal calibration checks. The RV templates and standards were 


reduced in the same way as the science target. 


^IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatories, 
which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative 
agreement with the National Science Foundation. 

‘ Muller et al. 1 2013ll show that the barycentric correction computed by the FEROS-DRS is only accurate to -^100 


m/s because it does not take into account coordinate precession. This small systematic error is not significant for 
our analysis, which achieves an RV precision of -^1 km/s. 
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2.3.8. CTIO-SMARTS 1.5-meter/CHIRON 


We obtained high-resolution (i?Ri28000) optical spectra of 2MASS J02155892-0929121 A and 
2MASS J06475229-2523304 A on the nights of UT 201 4 December 08 and U T 2014 December 
10, respectively, with the CHIRON echelle spectrograph (jTokovinin et alJl2013l i at the CTIO 1.5- 
meter telescope operated by the SMARTS Consortium at Cerro Tololo, Chile. CHIRON is an 
extremely stable fiber-fed instrument built for high precision (~1 m/s) RV exoplanet searches. Each 
spectr um covers the 4500 -8900 A spanning 62 orders. All observations are carried out in queue 
mode (jBrewer et al.ll2014l ) and are reduced in an automated fashion that includes bias subtraction, 
flat fielding, optimal spectral extraction, and wavelength calibration of each order using a ThAr 
exposure acquired after each target. 

We measure RVs for both targets by cross-correlating the science spectra with RV standards 
from IChubak et al.l (j2012l l that have similar spectral types and were also observed using the same 
setup with CHIRON. Each order is first normalized and the science and standard spectra are then 
linearly interpolated onto a common wavelength grid. After applying a barycentric velocity correc¬ 
tion to each spectrum, each order is cross correlated separately and the RV offset was determined 
by fitting a Gaussian to the cross correlation function. The final RVs (Table E]) are the mean from 
10 orders spanning 6120-6860 A, which were chosen to avoid telluric absorption features. Our 
adopted uncertainties incorporate both random errors (the standard deviation from the 10 orders) 
and an estimated systematic error of 0.5 km s“^ (measured by cross-correlating RV standards with 
each other) added in quadrature. 


2.3.9. Keck/ESI 


The wide separation of 2MASS J15594729-1-4403595 B from its host star (5"6) makes it 
amenable to spectroscopy without the need for AO corre ction. We observed 2MASS J15594729-I-4403595 
B using the Echellette Spectrograph and Imager (ESI; ISheinis et al.ll2002h on Keck H on UT 2012 
September 9 under near-photometric conditions. We used the UO slit, which yielded a spectral 
resolution of R ~ 4000 and wavelength coverage from 3000-10000 A. The slit was aligned perpen¬ 
dicular to the star-companion position angle to avoid contamination from the primary. We obtained 
three separate exposures, each 800 s, which yielded a S/N per pixel of ps 60 redward of 7500 A. 
Observations of the G191B2B and Feige 67 spectrophotometric standards were taken throughout 
the night. 


We reduced the data using the ESIRedux package (jProchaska et al.l l2003l : iBochanski et al 


2009 


fl 


ESIRedux performs bias, flat, and dark correction as well as extraction of the two- 
dimensional data into a one-dimensional spectrum. ESIRedux wavelength-calibrates the data using 


^http://www2. keck.hawaii.edu/inst/esi/ESIRedux/index.html 
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HgNe and CuAr arcs (obtained at the start of the night) to construct a two-dimensional wavelength 
map. The orders are combined into a single spectrum and an approximate flux calibration is applied 
(based on archived calibration data), which we improve upon using the observed spectrophotometr ic 
standards and a model of the atmosphere above Maunakea (|Buton et al.ll2012l : iMann et al.ll2013l i. 


3. Targets 

3.1. 2MASS 301225093-^2439505 AB 


_ 2MASS J01225093-2439505 B is a young mid-L dwarf companion identified bv iBowler et al 

((2013) during the PALMS high-contrast imaging planet search. Located at T.'5 (~52 AU) from the 


young active mid-M dwarf 2MASS J01225093-2439505 A, 2MASS J01225093-2439505 B possesses 
red near-IR colors and an angular iL-band spectral morphology, hallmarks of a dusty photosphere 
and low surface gravity, respectively. Our low-resolution optical spectrum of the primary from the 
SOAR/Goodman Spectrograph is shown in Figured! We find a spectral type of M4.0 ± 0.5 from 
visual and index- based classification with the IDL-based classification packag e Hammer (Table [3 


Covey et a 


Malo et al 


2007 1. which is slightly later than the M3.5V classification from iRiaz et al.l (120061 ). 
((2013) found the system probably belongs to the ~120 Myr AB Dor moving group, 
which was bolstered by Bowler et al. based on an RV measurement (9.6 ± 0.7 km s“^) and the 
photometric distance (36 ±4 pc). At this age, the companion has a luminosity that corresponds 
to both ~13 and ~25 Mjup masses based on substellar evolutionary m odels, which overlap in this 
region b ecause of delayed o nset of deuterium burning at lower masses ( Bowler et al.ll2013l L More 


recently. iGagne et al.l (j2014l l suggested t his system may be long to the 15-25 Myr j3 Pic moving group 
using the RV of 11.4 ± 0.2 km s“^ from iMalo et al.l (2013), in which case the companion mass is as 
low as Ri6 Mjup. A parallactic distance is needed to unambiguously distinguish these cases. Note 
that the RVs from Bowler et al. and Malo et al. differ by 2.5 a, indicating the p rimary could be a 


singl e-line spectroscopic binary (SBl). The Radial Velocity Experiment fRAVE: iKordopatis et al 
20131 ) found a slightly larger higher RV of 12.3 ± 2.5 km s“^, which is formally consistent with the 
other two measurements. 

On UT 2013 August 17 we obtained the first V-band images of 2MASS J01225093-2439505 B 
with NIRG2 (Figure [2|). Our astrometry is consistent with measurement s between 201 2 -2013 from 


Bowler et al. Based on the typical Y-J color of an M4.0V star from iRavner et al 


jiS), 


our 


V-band contrast of 7.6 ± 0.3 mag implies an apparent magnitude of V = 18.2 ± 0.3 mag (Table [3]), 
an absolute magnitude of My=15.4 ± 0.3 mag, and a Y-J color of 1.5 ± 0.3 mag. The relationship 
between My and spectral type has no t yet been mapped fo r young brown dwarfs and giant planets, 
but compared to old held objects from iDupuv fc Lid (120121 ). our new measurement best corresponds 
to L5-L7 objects. This is in good agr e emen t with our classihcation of L4-L6 based on H- and K- 
band spectroscopy from iBowler et al.l (j2013l i. 
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3.2. 2MASS J02155892-0929121 AabBC 


Riaz et al.l (200a) classify 2MASS J02155892-0929121 an M2.5 star displaying Ha emission 
{EW=-6.9 A) in their catalog of i?05AT-detected M dwarfs. We find a slightly later type of 
M3.5 with onr Good man spectrum (Table [Sj Figure [1]) but a similar Ha line strength of -5.6 A. 
Bergfors et al.l (120101 ) resolved this object into three components in a single epoch of imaging in 
2008 as part of the AstraLux Lucky Imaging survey: a close s tellar compan i on at C/.'6 (“B”) and a 
fainter component at 3'.'5 with red optical colors (“C”), which iJanson et al.l (|2012l l estimated has a 
spectral type of M8. 


Based on its activity and proper motion, iMalo et ahl (l2013lj suggest this sy stem is a possible 
member of the Columba, (3 Pic, or Tuc-Hor moving groups. iKraus et al.1 (j2014bl i also identify it as 
a candidate Tuc-Hor member, though their measured RV (10.1 ± 0.6 km s“^) disagrees with the 
ex pected value assum ing group membership by 5.2 km s“^. Three additional epochs were presented 


bv iMalo et ahl (120141 ): 0.5 ± 0.3 km s ^ and -0.6 ± 0.3 km s ^ in 2010, and 8.3 ± 0.3 km s 


..-1 


m 


2012. A consistent but less precise RV of 0.6 ± 5 km s ^ was measured by the RAVE survey in 
2007 (jKordopatis et al.ll2013l ). These RV variations indicate the primary is an SBl (Table [6]). 


We measure RVs of 6.2 ± 0.4 km s“^ and 8.3 ± 0.6 km s“^ for 2MASS J02155892-0929121 
from spectra obtained in 2012 December and 2014 December with Keck/HIRES and SMARTS 1.5- 
m/CHIRON. Our Keck/HIRES measurement agre es well with the predicted RV of 4.5 ± 1.4 km s“^ 


for Tuc-Hor membership from iMalo et al.l (j201^ ). Based on this RV, the web-based Bayesian 


tool BANYAN II to compute YMG membership probabilitiej^ by Gagne et al. ( 2014 1 returns a 
99.5% likelihood that this system belongs to the Tuc-Hor moving group. The weighted mean RV 
(4.0 ± 0.15 km s“^) is also consistent with the predicted value for Tuc-Hor. However, the changing 
RVs spanning 2010 to 2014 mean that more spectroscopic monitoring will be needed to derive a 
systemic RV and unambiguously determine its group membership. 

We imaged 2MASS J02155892-0929121 three times between 2012-2013 with Subaru/IRGS 
and Keck/NIRC2 (Figure |4]). We confirm the B and G components are comoving wi th the primary 


(Figur e [5]) and detect orbital motion between our epochs and those in 2008 from iBergfors et al 


( 201fll ). For our August 2013 data with NIRG2, only the central 198x248 pix^ (l'.'97x2'.'47) region 
of the detector encompassing the AB components was read out so we do not have astrometry of 
the G component at that epoch. 


3.2.1. Resolved AO Imaging of 2MASS ,102155892-0929121 AabB 

From visual inspection of our NIRG2 images in August 2012 and August 2013 we found that 
the A component appears slightly elongated compared to B at both epochs, but in the opposite 


®http://www.astro.umontreal.ca/~gagne/banyanII.php 
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directions (SE in 2012 and NE in 2013). The separation between A and B is close enough (C/.'6) 
that isoplanatic effects should be small, so any differences in the PSF shape are probably caused 
by a marginally-resolved close-in companion. 

Using a C/.'dxC/.'d region surrounding 2MASS J02155892-0929121 B as a reference PSF, we 
jointly fit two identical PSF templates to the A component w i th the Nelder-Me ad (“Amoeba”) 
downhill-simplex minimization technique (jNelder &: MeadI 1 19651 : iPress et al.l 120071 1 . Seven param¬ 
eters were allowed to vary: an x-position, y-position, and multiplicative scale factor for each PSF 
as well as an overall additive amplitude offset. Each image was fit separately (20 in 2012 and 10 
in 2013) and the resulting coadded frames are shown in Figure [6l In each case, subtracting the 
primary reveals a clear image of a companion and its first Airy ring. The joint fits provide excellent 
matches to the data as demonstrated by the low RMS residuals of ~8 DN (70 DN) compared to 
a peak flux of ~1600 DN (14000 DN) in our 2012 (2013) data. The median and standard devia¬ 
tion of the separations, position angles, and AT^-band flux ratios are 42 ± 7 mas, 112 ± 2°, and 
1.2 ± 0.3 mag for our 2012 data set. For our 2013 data we measure 42 ± 7 mas, 308 ± 5°, and 
1.17 ± 0.05 mag. 

Resolving the very close binary 2MASS J02155892-0929121 Aab at the same separation but 
with a nearly 180° change in P.A. suggests an orbital period of ~2 years (or shorter aliases). We 
perform a simple test to verify whether the amplitude of the R V variations are plaus ibly caused 
by this companion. Based on the empirically-derived relations in iBovaiian et al.l (|2012l ). the M3.5 
spectral type of the Aa component corresponds to an effective temperature of Ri3300 K. At an age 
of 40 Myr, this implies a mass of ~0.2 Mq for the primary (Aa) and ~0.09 for the companio n 
(Ab) using on the measured AT-band flux ratio and the evolutionary models of lBaraffe et al.l (Il998l l. 
Assuming a semi-major axis of ~1 AU, which roughly corresponds to the observed separation of 
the pair, a circular orbit, and an inclination of 90°, the induced velocity semi-amplitude is ~5.0 
km s“^. This is in good agreement with the maximum total RV amplitude change of ~11 km s“^ 
from random RV sampling over the past decade. 

There is no parallactic distance for this system, but spectrophotometri c distance estimates 
to 2MASS J021 55892-0929121 C (see below) and the unresolved primary ( Kraus et al. 2014b : 
Riaz et al.ll2006ll imply distances of fti20-30 pc. These correspond to physical separations of AU 
for Aab. Fortunately, the nearby components 2MASS J02155892-0929121 B and C (see below) offer 
excellent opportunities to measure individual dynamical masses of pre-main sequence M dwarfs 
using absolute astrometry of the system as opposed to a total mass with relative astrometry. 


3.2.2. The Ultracool Companion 2MASS J02155892-0929121 C 

Our 0.8-2.5 p.m. spectrum of 2MASS J02155892-0929121 C is shown in Figure [7l Overall the 
NIR SED agrees fairly well with the M7 field template. Individual bandpasses are most similar 
to M7-M8 objects. The H 2 O, H 2 O-I, and H 2 O -2 gravity-insensitive index-based classifications 
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from 


Allers &: Lid (j2013|) , which were originally defined by lAllers et alJ (j2007l l and ISlesnick et al 


(200J), give spectral types of M6.6 ± 1.7, M7.0 ± 2.3, and M6.7 ± 2.2 with a weighted mean 
of M6.7 ± 1.2. Altogether we adopt a final NIR spectral type of M7 ± 1. Of note are the 
shallow 1.244 ^m and 1.252 //m K 1 and 1.2 ^m FeH features in J band, hallmark sign atures 
of low surface gravity and youth (<200 Myr) in brown dwar f s (e.g ., iMcGovern et al.l 120041 ). We 


confirm this quantitatively using the indices from I Allers fc Liiil ( 20131 ). which yield an intermediate- 
gravity (INT-G) classification for 2MASS J02155892-0929121 G (Table [8]). This is strong evidence 
independent of kinematics that the system is quite young— likely ~30-150 Myr— implying that 
2MASS J02155892-0929121 G is indeed substellar. 

The iF-band p hotometric dis t ance t o 2MASS 402155892-0929121 C is 29 ± 5 pc using the 
field relations from iDupuv fc Lid (l2012l ). Since the system is young, this dis tance is probably 
underestimated. Nevertheless, adopting the AT-band bolometric correction from Ihiu et al.l ( 2010a ) 
gives a luminosity o f log (L/L^)=-3.23 T 0.16 dex, implying a mass of 42 ± 15 Mjup based on 
cooling models from iBurrows et al.l ( 1997 ). A larger distance of 40 pc corresponds to a luminosity 
of -2.95 ± 0.13 dex and a mass of 57 ± 19 Mj^p. 


3.3. 2MASS J02594789 0913109 AB 


This late -M companio n to th e M4.0 star 2MASS J02594789-0913109 A ( Reid et al.ll2007l ) was 
identified by iJanson et al.l (120121 1 from a single epoch of imaging in 2008. Aside from its large 
proper motion (/xtot=0'.T yr“^), the primary is otherwise unremarkable. It was not detected by 
ROSAT or GALEX and there are no signs the system is young. We do not detect lithium in 
our high resolution Kec k/HIRES spectru m and Ha is in absorption (i71T=0.34 A), implying an 
old age of several Gyr ( West et al.l 120081 ). Based on our RV measurement of 4.9 ± 0.5 km s“^, 
2MASS J02594789-0913109 is not consistent with any known YMGs at any distance. Our shallow 
imaging of the system with Keck demonstrates the pair’s shared proper moti on. At its est i mated 
distance of 36 ± 13 pc, an angular separation of Cf.'6 corresponds to ~22 AU. iJanson et al.l (j2012l l 
estimate a spectral type of M9.5 and a companion mass of ^ 0.08 We find a similar value 

of 87 ± 13 Mjup using a AT-band bolometric correction from iLiu et al.l (|2010al l. a uniform age 


distribution from 4-10 Gyr, and the lBurrows et al.l (|l997l l evolutionary models 


3.4. IRXS J034231.8-hl21622 AB (2MASS J03423180-hl216225 AB) 


IRXS J034231.8-1-121 622 is a young (60-30 0 Mvr: IShkolnik et al.ll2009l : IShkolnik et al.ll2012l ) 
M4.0 star at 23.9 T 1.1 pc (IDittmann et al.ll2014l l with Ha emission and saturated X-ray luminosity 
(jRiaz et al.ll2006l ). I.lanson et al.l (j2012l ) found a nearby (Cf.'S) candidate companion wi th red optical 


colors, which was confirmed to be a comoving brown dwarf and classified as LO ± 1 bv iBowler et al 


(120151 1 with near-infrared AO imaging and spectroscopy. Based on the age and luminosity of the 
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system, Bowler et al. inferred a mass of 35 ± 8 Mjup for IRXS J034231.8+121622 B. 

Here we present additional astrometry of IRXS J034231.8+121622 AB between 2007 and 
2013 with Keck/NIRC2 and Gemini-S/NICI. Combined with published astrometry (Table [Tj) there 
is clear orbital motion in the system. Although the coverage is quite limited, only a subset of 
Keplerian orbits are consistent with the observed astrometry. Our first attempt to constrain the 
orbital parameters and total mass of the system with orbit fitting is described below. 


The orbital elements are calculated by fitting the resolved astrometric observations to a Kep¬ 
lerian orbit using emcee (IForeman-Mackev et a,l 1120131 '). an implementation of the affine-invariant 
MCMC ensemble sampler of Goodman fc Wearel ( 201(1 ) . Our model has 8 free parameters: ^/ecosuJ, -y/e sino;, 
the time of periapsis, log(period), log(total mass), inclination, the position angle of the ascend- 
ing node, and distanc e. For the distance, we apply a prior following the parallax measured by 
Dittmann et al.l (|2014l h At each step, we forward model the orbits of the two components, calcu¬ 
lating their relative projected positions on the sky as viewed from the Earth at the time of each 
astrometric observation and determining the likelihood of these orbital elements (Figure [8|). We 
find that even though observations cover only 1.5% of the orbit, we are able to constrain the orbital 
period to 12% and the total mass to 18%. 

The median and 68% confidence intervals of the posteriors are 411 ± 51 years for the orbital 
period, 0.16 ± 0.07 for the eccentricity, and 0.152 ± 0.027 Mq for the total mass (Figure [9]). 

This is somewhat lower th an the total mass predic ted by evolutionary models. Using an 77-band 
bolometric correction from lCasaerande et al.l (|2008l ~) and the parallactic distance to the system, we 
find a luminosity of log L/Lq = -2.17 ± 0.05 dex for the primary star. This corresponds to masses 
of 0.15, 0.20, 0- 2 3, and 0.25 for ages of 60, 100, 200, and 300 Myr using the evolutionary models 
of lBaraffe et al.l (|l998l b The only way the total mass is consistent with these models is for young 
ages (<100 Myr), th ough it is also possible that systematic errors exist in the evolutionary models 
(jPupuy et al.l ()2014l )). Unfortunately, the relatively uncertain age of this system makes it difficult 
to use as a robust test of models. 


3.5. HD 23514 AB (2MASS J03463839+2255112) 


of 


HD 23514 is an F6 member of the Pleiades with an unusually high fractional infrared luminosity 
a2%, indicati n g the presence o f a lar ge amount of collisionally-ground dust within ~1 AU 


(jRhee et al.ll2008lb iRodriguez et al.l (j2012l l discovered a substellar companion to HD 23514 located 
at 2^6 (~360 AU) with AO imaging. HD 23514 B was also independently found bv I Yamamoto et al 


([2013) as part of the Subaru SEEDS exoplanet imaging program. As a companion to a well- 
established member of the Pleiades cluster, HD 23514 B is an excellent intermediate-age (120 ± 10 
Myr) benchmark for comparative analysis with field objects of unknown ages. 


Eigure [10] shows our H- and K- band (1.5-2.4 /rm) spectra of HD 23514 B taken with 
Keck/OSIRIS. We have not made any correction for reddening since extinction to the Pleiades 
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is small {E{B — V) < 0.06 mag; iTavlorl I 2 OO 8 I) . HD 23514 B exhibits typical features of late-M 
dwarfs including FeH absorption, strong CO bands, and Na 1 absorption. Interestingly, the H- 
band shape is only slightly angular, implying that late-M dwarfs showing even mor e pronounced 
angul ar morphology are younger than 120 Myr. This agrees with the conclusions bv lAllers &: Liu 
(|2013l b who found that intermediate- and very low-gravity brown dwarfs have ages <200 Myr. 
Compared to field templates, HD 23514 B best resembles M9-L0 objects in H band and M7--M8 
objects in K band. We therefore adopt a NIR spectral type of M8 ±1. In the future, J-band 
spectroscopy of HD 23514 B will be useful to better calibrate the age-dependency of J-band alkali 
absorption line strengths. 


3.6. G 160-54 AabBC (2MASS J04134585 0509049 AabBC) 


G 160-54 is a little-studied system with an (unresolved) spectral type of M4.5 ± 0.5 (Figured] 
Table [5|) and weak Ha emission {EW=-0.8 A from UH 2.2- m/SNIFS and EW=-0.bA A from 
CFHT/ESPaDOnS). Although not detected by ROSAT, a UV counterpart appears in GALEX 
{EUV = 21.7 =fc 0.4 mag. NU V = 20.79 ± 0.17 mag) and was resolved into a close visual triple 
system bv iBowler et al.l (120151 ) with NIR AO imaging at Keck and Subaru. G 160-54 B is a stellar 
companion 0 l'.2 (~4 AU) from the primary G 160-54 A and G 160-54 G is a fainter {AK=2.8 mag) 
comoving tertiary located at 3^.^3 (~70 AU. No parallax or RV has been measured for this system, 
but we note that at a distance of Ril8.5 pc and RV of «13 km s“^ its UVW kinematics ({-9.9, 
-20.7, -0.9} km s“^) would precisely match known members of the Tuc-Hor moving group. At that 
age (~30-40 Myr), the outer companion G 160-54 C would fall below the hydrogen burning limit. 

The very weak Ha emission, however, disagrees with established Tuc-Hor me mbers, which 


without exception have line strengths larger than -6 A at a spectral type of M4 ( Kraus et al 


2014bl ). Cross-correlating our high resolution spectrum of G 160-54 AB with RV standards yields 
three strong peaks, implying there is a fourth component in the system since G 160-54 C is too 
distant and too faint in the optical to produce such a large cross-correlation peak. Since we 
only have a single epoch we cannot identify whether this fourth member orbits G 160-54 A or B, 
but for this work we assume it orbits the primary “A” component in a similar configuration as 
2MASS J02155892-0929121. 

Our SpeX SXD spectrum of G 160-54 C is shown in Figure dH Gompared to field templates, 
it best resembles the M7 dwarf across the entire 0.8-2.4 fiui spectrum and among individual band- 
passes. There are no obvious signs of lo w surface gra v ity. T his is corroborated with index-based 
spectral and gravity classifications from lAllers &: Liul (|2013l i. which yield a NIR spectral type of 
M7 ± 0.5 and a field gravity (FLD-G). Altogether, despite the potentially coincidental kinematic 
overlap with Tuc-Hor, this system is certainly older than 150 Myr and most likely much older 
(>1 Gyr). At these old ages, all four components in this system are low-ma ss stars. Based on its 
distance (21 ± 9 pc), age (1-10 Gyr), and As-band magnitude (12.24 ± 0.06; Bowler et al. (2015)), 
G 160-54 G has a mass of 85 ± 16 using a A-band bolometric correction from Ibiu et al. ( 2010a ) 
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and the iBurrows et alJ (jl997l ) evolutionary models. The probability it is substellar (<75 Mjup) is 
only 7%. 


3.7. 2MASS J05464932 0757427 AB 


The M3.0 primary star 2MASS J054649 32-0757427 A is no t particularly noteworthy as it 
was not detected by ROSAT nor GALEX. iJanson et ahl (|2012l i identified a moderate-contrast 
companion (A2:'=5.6 mag) at 2'.'8 in two epochs of imaging in 2008 and 2009. We imaged the 
system in Ks band with NIRC2 in 2013. The known secondary was clearly detected and there are 
no signs of additional companions in the system. Janson et al. estimate a spectral type of >L0 for 
the companion based on z'- and i'-band contrasts. Our As-band contrast of 3.9 mag corresponds 
to an apparent magnitude of 13.7 mag and an absolute A^-band magnitude of RilO.4 mag at the 
estimated distance of ~45 pc from Janson et al. Assuming a 20% uncertainty in distance, the 
implied mass of 2MASS J0 5464932-0 7 57427 B is 87 ± 3 Mjup for an age ran ge of 1-10 Gyr using a 
bolometric correction from iLiu et al.1 (l2010al i and evolutionary models from IBurrows et al.1 (j 19971 1. 


3.8. 2MASS J06475229 2523304 AabB 


2MASS J06475229-2523304 is an active late-K-type star first recognized bv iRiaz et al.1 ( 20061 ). 
who found weak Ha emission {EW=-1A A) and a saturated X-ray fraction al luminosity (log Ly/Lhni 
3.13 dex). There are various conflicting classifications in the literature. iRiaz et al.1 (120061 1 find a 
spectral type of K7V from low-resolution optical spectroscopy . [Torres et al.1 (j2006l l find Klllle from 
high-resolution optical spectroscopy, and Pickles fc Deoagne ( 201ol l infer K3III from photometric 
fits to synthetic spectra spanning optical to near-infrared wavelengths. To better characterize this 
star we obtained a UH 2.2-m/SNIFS low-resolution optical spectrum spanning ~3500-8500 A. 


Figure [T^ compares our UH 2 .2-m/SNIFS spectrum to late-type templates from iBochanski et al 


( 20071 ) and I Mann et al.l ( 20131 ). 2MASS J06475229-2523304 broa dly resembles K7/M 0 dwarfs, and 
we find a consistent index-based spectral type of K7 using Hammer (ICovev et al.ll2nn7l ). We therefore 
adopt a K7.0 ± 0.5 classification for 2MASS J06475229-2523304. 

The band depths in our UH 2.2-m/SNIFS spect rum do not match dwarf templates in Figure [T^ 
however, which is confirmed via spectral indices from iReid et al.1 (jl995l l. We measure values of 0.890, 
0.950, and 0.952 for the Ti05, CaH2, and CaH3 indices, which sits above the dwarf locus in CaH2 
-|- CaH3 versus Ti05 diagram in a region occupied by giant stars (e.g., iMann et al.ll2012l ). This 
points to a low surface gravity, indicating 2MASS J06475229-2523304 is most likely a very young or 
evolved (subgiant or giant) star. The near-infrared colors of 2MASS J06475229-2523304 {J-H = 
0.59 ± 0.04 mag, H-Kg = 0.21 ± 0.04 mag) are similar to both late-K/early-M dwarfs and giants 
and cannot be used to distinguish between these evolutionary stages; only at later spectra. 1 types 
do their branches diverge in near-infrared color-color diagrams (e.g., iBessell &: Brettlll988l ). The 
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reduced proper motion of 2MASS J06475229-2523304 (ffy = V + 51og(/r) + 5, where /r is the total 
proper motion in arcseconds per year) is 10.3 mag, which is more consistent with field dwarf-like 
kinematic and photometri c properties instead of distant giant stars with high luminosities and low 
tangential velocities le.g.. iLeoine fc Gaidosll201in . 


Based on its activity and proper motion. iMalo et al.l ( 20131 ) suggest 2MASS J06475229-2523304 


compatible with the measured RV of -56.5 km s ^ from 

Torres et al. 

(2006). 

Malo et al. 

(2014) 

find an RV of -24.4 ± 0.8 km s ^ measured in 2011 and the RAVE survey ( 

Kordooatis et al. 

2013) 


find an RV of -67 ± 6 km s“^ measured in 2009, suggesting the primary is an SBl. Indeed, Torres 
et al. also listed this target as a possible binary (“SBl?”). We confirm this with our own RVs of 
-2.1 ± 0.9 km s“^ and -5.5 ± 0.8 km s“^ obtained in 2013 and 2014. Our du Pont spectrum from 
2013 shows slight Ha emission {EW=-ll.b A) with an asymmetric line profile; however, there are 
no signs of a companion in the cross-correllation function. Hereinafter this spectroscopic binary 
(SB) is referred to as 2MASS J06475229-2523304 Aab. 


There is no parallax measurement for this system vet. iRiaz et al.l (j2006l l find a J-band 


spectrophotometric distance of 45 ± 17 pc and I Ammons et al. (l2006l i estimate a distance of 
94^35 pc from SED fits. Assuming a Pleiades-like age, we infer an absolute R-ba nd magnitude 
of 8.1 1 ± 0.09 mag and a photometric distance of 36 ± 2 pc using the relation from lPowler et ^ 
( 2013 1. If 2MASS J06475229-2523304 Aab is younger than ~120 Myr or if it is an evolved star then 
its distance will be farther. Although we do not have a systemic RV, we note that no combination 
of RV between -50 to -2 km s“^ and distance between 10-100 pc is consistent with the XYZ 
spatial positions of known moving groups. 


3.8.1. The Comoving Companion 2MASS J06475229-2523304 B 

In 2012 we resolved a close (l^.^l), modest-contrast {AK=5.5 mag) point source near 2MASS J06475229- 
2523304 with Subaru/IRCS AO imaging (Figure [13]). Follow-up astrometry with Keck/NIRC2 in 
2013 and 2014 confirm the candidate is comoving with the primary star (Figure [14]). No orbital 
motion is detected. 

cite 60 mA Li I absorption from their high resolution spectrum of the 
primary star. Regardless of its kinematics, this would imply an age range of ~40-150 Myr assuming 
the system is indeed young rather than evolved, which raises the possibility that 2MASS J06475229- 
2523304 B is a new young brown dwarf companion. However, we do not detect any Li I absorption 
with our own high resolution du Pont/Echelle and SMARTS 1.5-m/CHIRON spectra, pointing to 
an age >100 Myr and a companion mass of about 40, 85, and 100 Mjup for distances of 40, 70, and 
100 pc, respectively. 

In December 2014 we obtained a iL-band spectrum of 2MASS J06475229-2523304 B with 
Keck/OSIRIS NGSAO. The spectrum does not resemble a late-M/early-L spectral type as expected 


Torres et al.l (2006j) 
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for a young (<200 Myr), nearby (within ~100 pc) system. Instead, the best match is an M3 
template with no obvious signs of low surface gravity (Figure [15]). This has forced us to re-evaluate 
the nature of this triple system since the component spectral types of K7 and M3 conflict with 
our measured NIR flux ratio of 5.5 mag. The corres ponding effective temperatures of K7 and M3 
dwarfs are 3960 K and 3412 K (|Bovaiian et al.ll2012l l. respect ively, which transla tes into a ilT-band 
contrast ratio of ~1.7-1.9 mag regardless of the system age (jBaraffe et al.lll998l ). In other words, 
2MASS J06475229-2523304 B is at least 3.5 mag fainter than it should be assuming all stars in the 
system are on the main or pre-main sequence. 


3.8.2. A Subgiant and M Dwarf Binary? 


Since it appears unlikely that 2MASS J06475229-2523304 AabB is composed of normal stars 
on the main or pre-main sequence, here we discuss a few alternative explanations. One option is 
that the system is very young and an edge-on protoplantary disk surrounding 2MAS S J06475229- 
2523 304 B is producing ^3.5 -4 mag of extinction in the NIR, simila r to TWA 30 B ( Looper et al 


20101 ) and possibly FW Tau (IKraus et al.ll2014al : iBowler et al.ll2014l ). However, this would require 
the system to be younger than the disk dissipation timescale of ~10 Myr, in which case we would 
expect lithium absor ption in the primary to be similar to TWA members (at least 300 mA; e.g., 
Mentuch et al.ll2008l ). This age is definitively ruled out with the nondetection of Li 1 in our high 
resolution optical spectrum, so this edge-on disk scenario is unlikely. 

Another option is that the primary is an evolved star while the compa nion is on the main 
sequen ce. Using the J-band contrast, the typical V-J color of M3 dwarfs fromiPecaut fc Mamaiek 
(j 20131 ) , and the absolute V -band magnitude of M3 dwarfs from [Drilling fc Landoltl (j2000|) , we find 
a spectrophotometric distance of 240 ± 40 pc to 2MASS J06475229-2523304 B. If we assume the 
primary is a red giant, its apparent U-band magnitude implies a distance of 1700 ± 100 pc. In¬ 
stead, to match the distance estimate of 240 pc for the companion, 2MASS J06475229-2523304 A 
would need to have an absolute U-band magnitude of ~4.0 mag. This is above the main sequence 
even if the pri mary is an equal- flux b inary and is m o st con sistent with evolved stars on the sub¬ 
giant branch. jPoimanskii (jl997l l and [Richards et al.l (j2012[ l measure a rotation rate of 7.56 days 
for 2MASS J06475229-2523 304 A which is somewhat faster than the typical rotation pe riod for 

How- 


subgiants of ~ 10-100 days (jdo Nascimento et al 


2012 


van 


Saders fc Pinsonneault 


201 


i). 


ever, the fast rotation p eriod can be explained by a tidally-locked short-period companion (e.g., 


Rvan fc Delivannis 19951 ). which can also account for its RV variations, Ha, and X-ray activity. 


We conclude that this system is most likely a distant (~240 pc) triple composed of an active K7 
subgiant primary with a close, likely tidally-locked companion and a much wider (~260 AU) M3 
dwarf. 





















































- 23 - 


3.9. 2MASS J08540240-3051366 AB 


This near by active M4.0 s tar was detected in ROSAT and GALEX and shows fairly strong 
Ha emission ( Riaz et al. 2006 . EW=-%A A). Its GALEX photometry {FUV=21A ± 0.4 mag, 
NUV = 19.53 ± 0.15) points to a ctivity levels compar a ble to known moving group members in 
UV/NIR color-color diagrams (e.g.. iBowler et al.ll2012bl ). IChauvin et al.l (l2010li estirnate an age of 


Eiaz et al. 2006). With 


100 Myr based on its saturated X-ray luminosity (log Ljy/Tboi=~2.79 dex; 
no RV, the BANYAN tool for calculating moving group membership probabilities ()Malo et al 


2 OI 3 I : iGagne et al.ll2014l ) predicts a high a priori value (87%) that 2MASS J08540240-3051366 is 
a member of /3 Pic. We measure an RV of 44.5 ± 0.6 km s“^ from our high resolution optical 
spectrum obtained in 2009. Based on this value, 2MASS J08540240~3051366 does not appear to 
be a member of a known moving group. 


Chauvin et al.l (j2010l ) first imaged this star as part of their AO imaging search for planets 
around young southern stars with VLT/NaCo. They identified a modest-contrast {AK=2.8 mag) 
companion located at l'!7. Assuming an age less than 200 Myr and a distance of ~ 10-15 pc as 
inferred by Riaz et al. and Ch auvi n et al.. the rnass o f 2MASS J08540240-3051366 B is below 
~50 Mjup. iJanson et al.l (j2012l ) and iJanson et al.l (j2014l ) present additional Lucky imaging of the 
system between 2010 and 2012 and find hints of orbital motion, which we confirm with our own 
NIRC2 images in 2012. 

In December 2014 we obtained Keck/OSIRIS H- and A-band spectroscopy of 2MASS J08540240- 
3051366 B. Compared to field templates in Figure fTOl it best resembles M4-M6 spectral types with 
no significant deviations from field objects that would point to an especially low surface gravity. 

It is therefore inconsistent with (3 Pic membership. We adopt a near-infrared spectral type of 
M5 ± 0.5, which is slightly earlier than the classification of ~M7 ± 2 suggested by our photometry 
in Table [3l 

Assuming an age of 100-200 Myr and a distance of 13 ± 5 pc, the implied mass of 2MASS 


J08540 240-3051366 B is 42 ± 21 Mt, ip based on the A-band bolometric correction from iLiu et al 


( 2010al ) and evolutionary models from iBurrows et al.l ( 19971 ). On the other hand, if the system is 
older then its mass is closer to the hydrogen burning limit (e.g., 78 ± 19 Mj^p at 0.5-1 Gyr). 


3.10. 2MASS J11240434-b3808108 AB 


This M4.5 sta r is a likely member of th e ^500 4= 1 0 0 My r Urs a Major moving g roup (see 


Bowler et al.l l2015l for a detailed summary). IClose et al.l (|2003l ) and ICruz et al.l (j2003l ) first no¬ 


ticed the late-M companion 2MASS J11240434-1-3808108 B located at 8"3 ( ^170 AU) from th e 
primary star. No additional companions were identified in deep AO imaging (jBowler et al.ll2015l ). 
We obtained a moderate-resolution spectrum of 2MASS J11240434-1-3808108 B with IRTF/SpeX 
(Figure [T7|l and optical AO imaging with Robo-AO (Figure fTHl) . We find a NIR spectral type 
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of M9.5 ± 0.5 using the spectral indices from lAllers &: Lid (j2013l l and vi sual compar i son t o field 
templates, which is somewhat later than the optical type of M8.5 found bv ICruz et al.l (|2003l i. The 
gravity indices of lAllers &: Liull2013l return a field surface gravity (FLD-G) with no signs of youth. 
This independently constrains the age of the system to >100 Myr, consistent with the tentative 
association to the Ursa Major moving group. 


3.11. PYC J11519+0731 AabB (2MASS J11515681+0731262 AabB) 

3.11.1. Spectral Classification and Binarity of PYC J11519+0731 Aab 


PYC J11519+0731 has received little previous attention i n the literature aside from id enti- 


hcation as a ROSAT source counterpart ( 

Zickeraf et al. 

2003: 

Haakonsen & Rutledee 

20091 and 

measurement of a 2.391 day rotationally-modulated photometric period ( 

Kiraeal 

20121. It was iden- 

tified as a candidate member of the fi Pic moving group by 

Schlieder et al. 

(2012b 1 from its proper 


motion and strong activity traced by ROSAT X-ray and GALEX UV emission. 


Schlieder et al.l ( 2012al l and lLeoine et al.l ( 20131 1 classify PYC J11519+0731 as M3 and M2.5, 


respectively, from moderate-resolution optical spectroscopy. Our Mayall/RC-Spectrograph optical 
spectrum and IRTF/SpeX near-infrared spectrum are shown in Figures [19] and [20l Even at the rel¬ 
atively low resolving power of our Mayall spectrum {R 2500), the Ha emission line {EW=-2.S A) 
is marginally resolved into two components (inset). This is easily split into a double-lined spec- 
troscopic binary at hig her resolutions (see below). Visual comparison to M dwarf templates from 


Bochanski et al.l (j2007l l indicates an M2 optical spectral type. Our spectrum deviates somewhat 
from the templates below ~7100 A, which is probably caused by slit loss from differential chro¬ 
matic refraction and/or the unresolved binarity. Altogether we a dopt a spectra l type of M2.0 ± 0.5. 
Compared to the NIR spectrum of the M2.0V star Cl 806 from Ravner et al. (2009) in Figure [20l 
PYC J11519-I-0731 shows some peculiarities in the form of deeper molecular and atomic absorption 
features primarily between 1.50-1.65 /rm, which could reflect its double-line spectroscopic binary 
(SB2) nature or a super-solar metallicity. Indeed, Schlieder et al. found rather large Na 1 A8200 
absorption (EW = 4.3 A), which is above the median value for field-age dwarfs of this type (see 
Figure 3 in ISchlieder et al.ll2012al ') and points to a high surface gravity and old age. 


We obtained three epochs of high-resolution optical spectroscopy of PYC J11519-I-0731 with 
Keck/HIRES, ESO-MPC 2.2-m/FEROS, and CFHT/ESPaDOnS between 2012 and 2014 (TableSj). 
Visual inspection of each spectrum reveals prominent emission and absorption features split into 
doublets (Figure El]), clearly indicating the primary star (PYC J11519-I-0731 A) is itself a double- 
lined SB caught with the components at a relatively large velocity (close physical) separation. 
Our HIRES spectrum exhibits strong, double-peaked Ha emission with dual self-absorption. The 
FEROS spectrum shows double-peaked Ca II H and K emission, strong Ha emission with dual 
self-absorption, dual reversed-emission cores in each of the Ca II IR triplet and the Na D lines, 
He I A5876 emission, and single-peaked H/3, Hy, and Hh emission. Similarly, the ESPaDOnS 
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spectrum displays double-peaked Ca II H and K emission, strong Ha emission with dual self¬ 
absorption, dual reversed-emission cores in each of the Ca II IR triplet and the Na D lines, dual 
He I A5876 emission, and dual-peaked H/3, H 7 , H(5 and He emission. 


We used used a cross-correlation (CC) package from iBender fc SimonI (120081 ) to determine the 
RV of the primary (Aa) and seco ndary (Ab) at each e poch, systemic velocity, and mass ratio of the 
system. Following the methods of iMazeh et al.l (j2002l 'l. we performed cross-correlations between the 
science spectrum and template spectra to measure the RV of each component. The templates can be 
observed standards or models, and the routine automatically accounts for the relative normalization 
of the spectra. 

The HIRES, FEROS, and ESPaDOnS spectra were analyzed in a similar fashion. The echelle 
orders were first trimmed to avoid telluric features, normalized, and, when necessary, the continua 
were flattened by fi tting and dividing by a polynomial. A set of K5-M4 RV templates were selected 
from iPrato et al.l (j2002l i and observed using EEROS. These same templates were used for the 
HIRES and ESPaDOnS analysis. We attempted CC with single, binary, and rotationally-broadened 
templates. Eor our FEROS data, the strongest CC power (~90%) was achieved for the primary and 
secondary using an M2.5-I-M2.5 system with i?I/p=15.7 ± 1 km s“^ and vpsinip = 14 ± 2 km s“^ 
and RVs = -18.2 ± 1 km s“^ and vssinis = 12 ± 2 km s“^. The RV uncertainties are dominated 
by the use of observed template spectra, and the projected rotational velocity uncertainties reflect 
the velocity increments used to generate the suite of broadened templates. For our HIRES data we 
measure CC peaks at RVp=50.9 ± 1 km s“^ and RVs = -51.9 ± 1 km s“^, and for our ESPaDOnS 
data we find RVp=69.5 ± 1 km s“^ and RVs = -75.5 ± 1 km s“^. The range of flux ratios we 
infer are between 0.5-0.75 at ~6400 A. 

With three epochs of RV data for each component, th e mass ratio of the system (q) and 
systemic RV ( 7 ) can be determined following the methods of IWilsonI (jl94ll ). Figure [22] shows the 
component velocities plotted against each other; the mass ratio of the system is the negative of 
the slope of the line. We find a nearly equal mass ratio of q = 0.95 ± 0.01 and systemic velocity 
7 = -0.3 ± 1.0 km s“^. 

The large velocity semi-amplitude of PYC J11519-I-0731 Aab implies a short orbital period 
(at most a few days), so the high ene rgy emission is most likely a result of fast rotation from 
tidal synchronizat ion rather than youth ([Torres et al.ll2002l : [Torres et al.ll2003l : [Shkolnik et al.[[20iril : 
Kraus et al.l[201l[ ). Further RV and photometric monitoring could provide a spectroscopic orbit 
and rotationally-modulated period, which would give a better indication of whether tidal spin up 
could be responsible for the observed activity. 


3.11.2. Age and Kinematics of PYC J11519+0731 Aab 

In general, Li I absorption at 6708 A from either component of PYC J11519-|-073lAab would 
provide strong, independent evidence of youth. The top panels of Figure!^ show the Li I region in 
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our three high-resolution spectra. Our HIRES spectrum shows an absorption line that is close to 
the predicted position of Li I in the primary. However, in the other epochs there is no convincing de¬ 
tection of Li I in either component. At the 12-25 My r age of the j3 Pic moving group (|Yee &: Jensen 
2010 : iBinks &: JeffriesI l2014l : iMamaiek h Belli l2014l i , lithium is not typically detec ted for spectral 


types of M2; only at younger ages of <10 Myr is photospheric lithium still present (jMentuch et al 
20081 ). The lack of strong lithium in either component sets a lower age limit of ~10 Myr for the 
system and is therefore consistent with the age of /? Pic. 

Despite lacking a parallax, our systemic RV enables constraints on the kinematics of PYC 
J11519-I-0731 Aab. Figure 1231 shows the range of UVW heliocentric space velocities (based on 
the procedur e from [Johnson &: SoderblomI 119871 with the coordinate transformation matrix from 
Murray! 1 19891 ) and XYZ galactic positio ns for PYC J11519-I-0731 compared to members of young 
moving groups from Torres et al. ( 2008 )^1 The systemic RV we measure (- 7 = -0 .3 ± 1.0 km s“^) 
agrees with the value of -1.5 ± 1.0 km s“^ predicted bv ISchlieder et ah! ( 2012bl i assuming /3 Pic 
membership, and the range of UVW velocities is consistent with that group for distances between 
~20-30 pc. However, the physical location of PYC J11519-I-0731 disagrees from known members 
in Z by at least 20 pc. 

Group membership can be further tested by placing PYC J11519-I-0731 on a color-magnitude 
diagram. Because they are still contracting, (3 Pic memb ers are over lumino us by ~1 mag in V 
band compared to field stars on the main sequence (e.g., iBiedel et al.ll2014l ). Assuming similar 
spectral types for PYC J11519-I-0731 Aa and Ab and a near-IR flux ratio of 0.6 ± 0.1 from our 
CC analysis ( Sectio n 13.11~T]) . the integrated apparent V-band magnitude of 12.5 ± 0.1 mag from 
Droege et ah! (j2006li can be decomposed into individual V-band magnitudes of 13.01 ± 0.12 mag 
and 13.57 ± 0.16 mag for Aa and Ab. At a distance of 20 pc (30 pc), these correspond to absolute 
magnitudes of 11.50 ± 0.12 mag (10.62 ± 0.12 mag) and 12.07 ± 0.16 mag (11.19 ± 0.16 mag) for 
Aa and Ab, respectivel y. At the system V -K color of 4.6 ±0.1 mag, (3 Pic members are expected 
to have My^^.b mag ( Riedel et al. 2014). Instead, both components of PYC J11519±0731 Aab 
appear consistent with the main sequence given the range of plausible distances assuming group 
membership. It is therefore likely this system is simply an old but active interloper with (partially- 
constrained) kinematics similar to the f3 Pic YMG. 


^Note that the pr oper motion of PYC J11519-I-0731 is incorrectly listed in several all-sky catalogs. USNO- 
B I Monet et al. 20031 ') and the updated PPMXL Catalog i Roeser et al. 201ol ) list proper motions of ^aCosS = - 
112 mas yr“^ and fis = 122 mas yr“^ and fj.aCosS = -126 mas yr“^ and fis = 110 mas yr“^, respectively. However, 
visual inspection of older Digitized Sky Survey and 2MASS images clearly shows movement in the southwest direction. 
We therefore adopt the original PPMXL proper motion of fj,aCos5 = -142 ± 19 mas yr“^ and jj.s = -89 ± 18 mas yr“^. 
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3.11.3. PYC J11519+0731 B 


In 2012 we imaged PYC J11519+0731 as part of our PALMS survey. As expected, the close SB2 
pair PYC J11519+0731 Aab was not resolved in our NIR images, implying a separation less than 
about one PSF FWHM (40.4 ± 0.8 mas in H band, or ftil.5 AU at 37 pc). Indeed, two 0.2 Mq stars 
on circular, edge-on orbits would need to reside at separations of 0.03 AU to match the maximum 
observed RV semi-amplitude. We did, however, identify a moderate-contrast (A77=5.4 ± 0.2 mag) 
point source located at psC/.'S from PYC J11519-1-0731 Aab and immediately obtained unsaturated 
frames in Y, J, H, Ks, and L' (Figure [241 and Table [3]). The companion’s red colors point to 
a cool temperature and late-M spectral type. To search for additional companions in the system 
we obtained deeper 77-band coronagraphic images in 2013 (20 frames each with 30-sec integration 
times) in ADI mode. 


A single image in our ADI sequence is shown alongside the PSF-subtracted and coadded version 
using the entire sequence in Figure!^ The companion is visible in the single image and is similar 
in appearance to the quasi-static speckles, but stands out prominently with a peak signal-to-noise 
ratio of 370 in our PSF- subtracted irnage. T he PSF subtraction was carried out in a similar fashion 
as the deep imaging in Bowler et al. (2 015 Il except h ere we have made use of 100 reference PSF 
images for our LOCI reduction ([Lafreniere et al.l 120071 1 from a PSF library of over 2000 registered 
coronagraphic images from our PALMS survey. Reference images were selected as having the lowest 
residuals between C/.'3-2" after scaling and subtracting the library template to each science image. 
Local masking was applied at the location of PYC J11519-I-0731 B to minimize biasing the LOCI 
coefficients. 


We obtained three additional epochs of astrometry with Keck/NIRC2 between 2012 and 2014, 
and one with Gemini-S/NICI in 2013 (Table [3]). Compared to the expected motion of a stationary 
object fFigure [26|) . our astrometry unambiguously shows that PYC J11519-I-0731 B is comoving 
with the primary Aab. We also find clear signatures of orbital motion. The reduced chi-squared 
iXu) value of a constant (linear) fit to the separation measurements is 7.7 (1.0), and a constant 
(linear) fit to the PA is 268 (3.5). The linear fits correspond to changes of -1-9.7 ±1.2 mas yr“^ 
and “2.04 ± 0.04 ° yr“^ for the separation and PA, respectively. Curvature in the orbit should be 
evident within a few years. 


Our AT-band Keck/OSIRIS spectrum of PYC J1 1519±0731 B is shown in Figure!?^ comp ared 


to field templates from the IRTF Spectral Library (ICushing et al.l l2005l : iRavner et al 


2009|) and 


low-gravity (VL-G) templates from lAllers &: Liul (j2013l L PYC J11519±0731 B appears most similar 
to the held M8 template in overall shape and the depth of CO and Na I features. The 2.29 fiui CO 
absorption is weaker and the ~2 /rm steam absorption is stronger in the very low gravity objects 
compared to PYC J11519±0731 B. In addition, the 2.207 fim Na I doublet, which is sensitive to 
surface gravity (and hence age), is weaker in all very low-gravity templates. This is perhaps the 
strongest indication that PYC J11519±0731 B is older than ~30 Myr. Because of the limited 
spectral coverage, we adopt a spectral type of M8 with uncertainty of ±1 subtype. 































- 28 - 


The mass of PYC J11519+0731 B depends on the distance and, more importantly, on the sys¬ 
tem age. The photometric distance to PYC J11 519-I-0731 B using the 77-band absolute magnitude- 
spectral type relations from lPuDuv &: Liul (j2012l l is 37 ± 6 pc. This incorporates uncertainties in the 
phot ometry, spectral c lassification, and rms of the relation. Using the Tf-band bolometric correction 
from iLiu et al.l (|2010al l we infer a luminosity of logL/L*=-3.35 ±0.16 dex. A t ages of {21 ± 2 Myr , 
100 ± 10 Myr, 1.0 ± 0.1 Gyr, 5.0 ± 0.5 Gyr}, the evolutionary models of Burrows et al. ( 1997 ) 
imply masses of {15.3 ± 1.3, 42 ± 7, 89 ± 5, 90 ± 4} Mjup for PYG J11519±0731 B. If the system 
is younger than ~350 Myr then the companion mass lies below the hydrogen-burning limit. 

No other candidate companions were identified in our deep imaging, which extends out to ~4'' 
for complete azimuthal sensitivity and to ~8^^ for partial sensitivity. The 7 a 77-band contrast curve 
shown in Figure [2^ is measured by computing the rms in 5 pixel (50 mas) annuli as described in 
Bowler et al. ( 20151 ). Unsaturated frames are used at separations below 0'.'3 and our deep ADI data 
are used at larger radii. Our hnal contrast curve in AH and held of view coverage (in brackets) 
is {4.0, 5.4, 10.6, 12.2, 13.1, 13.8, 13.4, 14.0, 14.2, 14.0, 13.4} mag for separations of {O'.'!, 0'.'2, 
0'.'5, 0'.'75, 1"0, 1"5, 2", 3", 4" [1.0], 5" [0.69], 8" [0.05]}. That is, we are fully sensitive to spatial 
coverage within ~4", but are only sensitive to a fraction of the sky beyond that owing to the square 
detector. At ages of {21 Myr, 100 Myr, 1 Gyr, 5 Gyr}, we remain highly sensitive (>90%) to 
co mpanions with mass es as low as {1, 2, 10, 25} Mjup at 100 AU based on the evolutionary models 
of Baraffe et al. (2003). 


3.12. G 180-11 AB (2MASS J15553178±3512028 AB) 

A T.T companion to the M4.0 ± 0.5 star G 180-11 (a.k.a. GJ 3928) was hret di s covered by 


McGarthv et al. ( 

20011 and confirmed to be physically bound by 

Daemeen et al. 

(20071 

. Law et al. 


( 20081 ), 


and 


.lanson et al.l (|2012l ). Unr esolved detectioiis in ROSAT and G ALEX together with 


strong Ha emission (PTUss-O t o -8.1 A; Riaz et al.ll2006l : IShkolnik et al.ll2009l ) indicate the primary 


and/or companion are active. Riaz et al. ( 2006l l found a saturated fraction al X-ray luminosity o f 


log Ljv/Lboi=“2.98 dex for the system, pointing to a relatively young age. iHartman et al.l (j201l|) 
measured a rotational period of 3.5209 days as part of the HATNet survey. 


Recently, iDittmann et al.l (120141 ) measured a parallactic distance of 28 ± 3 pc to G 180-11 AB 
as part of the MEarth transit survey. At this dista .nce. the primary-corn panion separation is 45 AU. 
Together with the RV of-15.5 ± 0.7 km s“^ from I Shkolnik et al. ( 2012 1. this corresponds to UVW 
kinematics of {-35 ± 3, -17.6 ± 1.4, and 8 ± 2} km s ~^ and XYZ space positions of {9.9 ± 1.1, 


20131 . 


14.9 ± 1.6, and 21.5 ± 2} pc. iMalo et al.l (j2013l . l2014l l hnd that G 180-11 AB is a likely member 


of the Argus moving group based on the RV alone, which matches the predicted value assuming 
group membership, but the parallactic distance disagrees with the statistical distance of 13 ± 1 pc 
so the similar RV is probably coincidental. Indeed, the kinematics of G 180-11 AB do not appear 
to be consistent with any known moving groups. 
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The pair are easily resolved with our Keck/NIRC2 and P60/Robo-AO obs ervations spanning 
0.6-2.2 /xm (Figure [18]) Compared with previous astrometry for this system Isee lJanson et al.ll2012 
for a summary), our new data show slight outward motion by ~6 mas yr“^ and P.A. change of 
~0?8 yr“^. The orbital motion is quite slow so it will take decades to complete an appreciable 
portion of its orbit. 


Daemgen et a 

1. (2007 

) infer 

). Cruz et al. ( 

2003 

) deduce 


ratio, I Cruz et al.l (|2003l i deduce an ~M6 classihcation, and iLaw et al.l (|2008ll estimate a spectral 
type of M6.5 from its optical colors. In 2013 we obtained resolved near-infrared (JHK) spec- 
tra of G 180-11 B with Keck/OSIRIS. Compared to templates from the IRTF Spectral Library 
( Ravner et al. 2009) in FigureG 180-11 B best resembles M5-M6 field objects from 1.2-2.4 /xm 
and in individually normalized J, H, and K bands with a slightly better match to the M6 spec¬ 
trum. We therefore adopt a near-infrared spectral type of M6.0 ± 0.5, in agreement with previous 
estimates from colors. Although the 2.21 /xm Na 1 doublet line is equally strong in G 180-11 B as 
in the field templates, the 1.244 and 1.253 /um K I lines are noticeably weaker, pointing to a low 
surface gravity in line with an age younger than ~200 Myr. 


Based on the parallactic distance and an Lf-band bolometric correction from iLiu et al.l (|2010al L 
we measure a luminosity of log L/Lq = -2.48 ± 0.13 dex for G 180-11 B. At ages of 50 ± 5 Myr, 
100 ± 10 Myr, and 200 zfc 20 Myr, this c o rresp onds to masses of 88 ± 18 Mjup, 120 ± 20 Mjup, and 


150 ± 20 Mjup using the lBurrows et al.l (jl997l ) evolutionary models. Unless this system is younger 
than ~40 Myr, G 180-11 B is a low-mass star instead of a brown dwarf. 


3.13. 2MASS J15594729-b4403595 AB 


2MASS J15594729-b4 403595 is a young ac tive early-M dwarf with a comoving substellar com¬ 
panion first identified by jjanson et al.l (|2012l . see I Bowler et al.l l2015l for details) at a projected 
separation of 5'.T (150 AU). The companion was easily recovered with our Robo-AO imaging in i' 
and z' bands f Figure fTHI) . Our SNIFS 4000-9000 A optical spectrum of the primary in Figured] 
reveals weak Ha emission (FiVF=-2.3 A) an d is well-matched with an M1.5V spectrum, created 
by averaging the Ml and M2 templates from iBochanski et al.l (j2007l ). Our o wn RV of -19.6 zb 0.6 
km s“^ from ESPaDOnS disagrees with the value of -15.8 ± 0.5 km s“^ from iMalo et al.l ( 20141 ) by 
4 km s“^, a 4.9 cr difference (TableE]), raising the possibility that the primary is an SBl. Although 
the system does not appear to coincide with any known moving groups, the intermediate-gravity 
near-infrared spect rum of 2MASS J15594729-I-4403595 B implies a young age of Ri50-200 Myr 
( Bowler et al. 2015 1. 


Our spatially-resolved Keck/ESI optical spectrum of 2MASS J15594729-I-4403595 B is shown 
in Figure (301 2MASS J15594729-I-4403595 B exhibits modest Ha emission {EW=-QA A) and 
absorption features typical of late-M dwarfs. Compared to templates from Bochanski et ah, 
2MASS J15594729-I-4403595 B fits the M7.0 and M7.5 objects reasonably well. The index-based 
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classification from Hammer (|Covev et al.ll2007l i is M7. Altogether we adopt an optical classification 
of M7.5 ± 0.5. Note, however, that the best match to the entire 0.5-2.5 ^m spectrum is slightly later 
(M8: Figure [3T]) . which is common with visual classification of young brown dwarfs (jAllers Liu 

2OI3I 1. 

As shown in Figure IHTl we find strong Li I A6708 absorption in the companion {EW = 0.71 A). 

The detection of lithium pr ovides independent constraints on the mass and age of the com p anion 
throug h the “lithium test” ( Rebolo et al. 1992 : Rebolo et al. 1996 : Kirkpatrick et al. 2008 1. Basril 
(| 19981 1 shows that an object with Li I absorption with a luminosity below log (L/Lq) = -3.05 dex or 
Tgff < 2670 K (spectral type >M7) is unambiguously substellar and young (<200 Myr). This is con¬ 
sistent w ith the the int e rmed iate surface gravity and mass of 43 ± 9 Mjup for 2MASS J15594729+4403595 B 
found in Bowler et al. ( 2015 1. 


3.14. GJ 4378 Aal2b (2MASS J23572056-1258487 Aal2b) 

3.14-1- The wide binary GJ 4378 A and GJ 4379 B 

GJ 4378 A (LP 704-15, G 273-186) and GJ 4379 B (2MASS J2357 1934-1258406, G 2 73-185, 
LP 704-14) form a pair of comoving mid-M dwarfs separated by 20" (iGiclas et al.lll975ll. Both 
stars have similar reported K-b and (13.0 mag from A PASS/UGAG4, IZacharias et al.ll2013l l and 
r'-band (12.4 mag from GMG15, iMuihos &: Evansll2014ll m agnitudes, but GJ 4378 A is brighter in 
JHK by ^0.5 ma g (from 2M ASS, Skrutskie et al. 2006 1. GJ 4378 A is classified by Reid et al 


(200J), iReid et al.l ( 1995l l. and iRiaz et al.l ( 2006l l as M4, M4, and M3, respectively, from optical 
spectroscopy. The same studies classify GJ 4379 B as M4, M3, M4. We find similar spectral types of 
M3.5 for GJ 4379 B and M4.5 for GJ 4378 A from our own low-resolution optical spectra (Figured]). 
Although GJ 4378 A is slightly brighter, its spectral type is slightly later. Our moderate-resolution 
SOAR spectra of the pair are shown in Figured GJ 4378 A exhibits Ha emission (FTT=-4.5 A) 
while none is apparent in the companion GJ 4379 B. 

The system coincides with a ROSAT Bright Source Gatalog detection (IRXS J125720.0- 
125852), but it is unclear which compone nt is the counterp art because of the ROSAT PSPG’s 
large positional uncertainty of ~13" (1 a; IVoges et al.lll999l l. GJ 4378 A was detected in both 
GALEX in FUV and NUV bandpasses, while GJ 4379 B was detected in NUV alone. The 
NUV-J color, which traces UV excess from chromospheric emission, is 10.58 ± 0.05 mag and 
13.17 ± 0.19 mag for the A and B, respectively. The blue NUV-J of GJ 4378 A is consistent with 
activity levels of YMG members (10-150 Myr), b ut the red color of GJ 4379 B is more consistent 
with older M dwarfs in the Hyades and the field (|Rodriguez et al.l 120131 ). 


Malo et al.l (j2013l ) identify GJ 4378 A as a candidate member of the ~40 Myr Argus asso¬ 
ciation based on its proper motion and activity. Assuming group membership, they predict a 
kinematic distance of 25 ± 1 pc and RV of -1.0 ± 0.8 km s“^. In 2009 we obtained echelle spec- 
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troscopy of GJ 4378 A and GJ 4379 B with the du Pont telescope. Cross-correlation analysis of 
our du Pont spectrum of GJ 4378 A shows it is a near-equal flux SB2 with component velocities 
of -7.0 ± 0.8 km s“^ for the primary (GJ 4378 Aal) and 49.2 ± 1.1 km s“^ for the companion 
(GJ 4378 Aa2jl|. Assuming the pair are equal mass, the systemic RV of GJ 4378 Aal2 is 21.1 ± 0.7 
km s“^. We measure an RV of 20.6 ± 0.4 km s“^ for GJ 4379 B and find no signs it is a spectro¬ 
scopic binary. Assuming the RV for GJ 4379 B represents the systemic velocity, the range of UVW 
values for distances of 20-60 pc do not coincide with any known moving groups. 


Despite their similar spectral types and presumed coevality, GJ 4378 A and GJ 4379 B appear 
to have very different chromospheric activity levels. Like PYC J11519-I-0731 Aab (Section 13.111) . 
GJ 4378 Aal2 is an SB2 with enhanced UV emission (and the likely source of the X-ray emission) 
probably induced by tidally locked stars instead of intrinsic fast rotation from youth. On the other 
hand, GJ 4379 B appears to be single and is therefore a better tool to diagnose the system age. 
Indeed, neither our spectroscopy nor shallow Keck /AO iL,c;-band ima ging of GJ 4379 B on UT 2014 
August 04 shows any signs of a close companion. IWest et al.l (j2008l ) find a characteristic activity 
lifetime of 4.515 'q Gyr for M4 stars as traced by the presence of Ha emission. The lack of Ha 
emission in GJ 4379 B therefore suggests an old age of >4 Gyr for the system. 


3 . 14 . 2 . GJ 4378 Ab 


We imaged GJ 4378 A in Oct 2012 with IRGS on Subaru, easily resolving it into close (0'.'47) 
modest-contrast {AH=3.1 mag) binary. Astrometric monitoring throughout 2013 and 2014 confirm 
the pair are comoving (Figure [3^ and undergoing significant orbital motion of C/.T564 ± C/.'OOOO 
yr“^ in separation and -1?06 ± 0?02 yr“^ in PA. 


In August 2013 we obtained 1-4 /rm Keck/NIRC2 AO images of the system (Figure 1331) . 
Decomposed photometry of the pair presented in Table [3] imply the companion has red colors 
consistent with a late-M dwarf. This is confirmed with our H- and K-band OSIRIS spectroscopy 
of GJ 4378 Ab (Figure [33]), which shows atomic and molecular absorption features c haracteristic 
of late-M/early-L dwarfs. Gompared to field late-M dwarfs from [Gushing et al.l ((2005[), the overall 
1.5-2.4 ^m spectrum is most similar to the M9V template. The 77-band spectrum alone provides 
a reasonable match to the M9V template, while the 77-band spectrum is an excellent match to the 
M8V template. Altogether we adopt a spectral type of M8 ± 1 for GJ 4378 Ab. Interestingly, the 
77-band spectral shape appears somewhat more angular than the field objects, which would imply 
a young age (<150 Myr) for the system ([Allers &: Lid [2013[ 1 and appears to contradict the old 
age inferred from the lack of Ha emission in the wide companion GJ 4379 B. Moderate-resolution 
J-band spectroscopy of GJ 4378 Ab would help address the system age since the 1.1-1.3 /rm region 


^To avoid confusion, we adopt the following nomenclature for this hierarchical quadruple system: we retain the 
names of the 20” pair GJ 4378 A and GJ 4378 B, GJ 4378 A is a visual binary fsee I3.1472ll with GJ 4378 Aa 
and GJ 4378 Ab components, and GJ 4378 Aa is an SB2 with components GJ 4378 Aal and GJ 4378 Aa2. 
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possesses a number of gravity-sensitive features. 


Using //-band spectral type- absolute rnagnitud e relations from IDuduv &: Liul (1201211 and H- 
band bolometric correction from iLiu et al.l (j2010al l. we infer a spectrophotometric distance of 
14 zb 2 pc for GJ 4378 Ab and a luminosity of logfLhni/ U*)=-3.35 zb0.14 dex. This agrees with 
the photometric distan ce of 12 pc to GJ 43 78 Aa from iLepine fc GaidosI (j2nill ). Based on the 
evolutionary models of I Burrows et al.l (|l997l l. the luminosity of GJ 4378 Ab implies a mass of 
90 zb 3 Mjup for an age of 5 zb 2 Gyr with little dependence of age beyond ssl Gyr. 


4. Discussion and Conclusions 

The high density and sheer number of M dwarfs in the solar neighborhood imply that many old 
stars will share similar UVW space motions with nearby YMGs. The task of sifting through this 
haystack of potential members is made even more difficult when only an RV or parallax is available, 
and in most cases neither of these are. We have undertaken a detailed kinematic and spectroscopic 
analysis of 13 K and M-type systems hosting low-mass companions to clarify their ages, companion 
masses, and possible associations with YMGs (Figure [35]). Eleven of the primary stars have been 
previously identified as being young in the literature, eight of which have been listed as candidate 
members of YMGs. Three of the late-type companions to moving group candidates were discovered 
during our PALMS AO imaging survey and are new in this work (2MASS J06475229-2523304 B, 
PYC J11519-b0731 B, and GJ 4378 Ab). 

Only two of the systems appear to be likely members of YMGs (2MASS J01225093-2439505 
AB and 2MASS J02155892-0929121 AabBC). Another four show signs of youth but are probably 
not kinematically associated with any known moving groups (IRXS J034231.84-121622 AB, 2MASS 
J08540240-3051366 AB, G 180-11 AB, and 2MASS J15594729-b4403595 AB). With the exception 
of 2MASS J11240434-b3808108 B, whose kinematics are consistent with the ~500 Myr Ursa Major 
moving group, the remaining systems with low-mass companions are probably old. In particu¬ 
lar, we have shown that the host stars to 2MASS J06475229-2523304 B, PYC J11519-b0731 B, 
and GJ 4378 Ab are close spectroscopic binaries with little or no evidence that they belong to 
nearby young moving groups. This underscores the need for follow-up RV (and ideally parallax) 
measurements for new proposed members of YMGs to identify old M-bM SBls and SB2s that can 
masquerade as young stars. 

2MASS J02155892-0929121 AabBC is likely to be an especially useful system in the future. If 
its systemic RV and parallax confirm it is a member of Tuc-Hor, which is suggested by the mean 
(4.0 zb 0.15 km s“^) and range (-0.6 to 10.0 km s“^) of RV measurements, 2MASS J02155892- 
0929121 will offer an excellent opportunity to measure dynamical masses for the Aab components 
and, eventually, the B component. Since the outer companion (“C”) is a young brown dwarf, this 
quadruple system offers the opportunity to test pre-main sequence evolutionary models across the 
substellar boundary. Regardless of any possible age spread in Tuc-Hor itself, the components of 
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this system were probably formed simultaneously, removing any ambiguity about pos sible age dif- 
ferences. This is analog ous to the well-stu died young high-order be nchmarks LkCa 3 (ITorres et al 


2013 1. HD 98800 (e.g., Prato et ahlbooil l. and especially GG Tau ( White et al. 1999 1. which also 
contains a subs tellar member and re cently gained a fifth component from sparse aperture masking 
interferometry (jDi Folco et al.ll2014^ . 


Four systems in our sample are examples of close spectroscopic binaries on what are likely 
tidally-circularized, few-day orbits with nearby lower-mass tertiaries that we have directly re¬ 
solved with AO (G 160-54 AabBG, 2MASS J06475229-2523304 AabB, PYC J11519-F0731 AabB, 
GJ 4378 Aal2b). Spectroscopic binaries do not form in situ but instead must lose angular mo- 


(e.g., 

Gorti & Bhatt 

1996: 

Sta 

ilei 

20161 or through dynamical interactions with a third body (e.g.. 

Eesfleton & Kiseleva-Eesrleton 

2001 

). These close hierarchical triples are consistent with this latter 


scenario in which the tertiary remained relatively close in after a dynamical encounter with the 
(now tid ally lockedl secondary star instead of being ejected to a wide orbit or completely out of th e 


system (ITokovinin et al.l 120061 : iBeipurth et al. 


2010 


Allen et al.l l2012l : iReipurth fc Mikkolal l2013l l . 


These tertiaries are expected to have relatively high eccentricities, which can eventually be tested 
with continued astrometric monitoring. 

Regardless of their ages, many of these systems are usef ul for calibra t ing em pirical metallicity 
relationships at the bottom of the main sequence. Recently iMann et al.l (120141 ) presented an em¬ 
pirical method to measure the metallicities of late-M dwarfs with A'-band spectra by leveraging the 
metallicity of wide, earlier-type binary companions for calibration. Although this was the largest 
such study to date, only a handful of late-M dwarfs in binaries were used (7 systems beyond M7.5), 
so building larger samples of these benchmark systems will help better constrain these relationships. 

Gaia will soon provide precise parallaxes for the currently unknown stellar members of YMGs, 
but its estimated 1-15 km s“^ RV precision is generally not adequate to differentiate moving 
groups with UVW space velocities on ly a few km s“^ apart like TWA, Tuc-Hor, and Carina 
( Lindegren et al.ll2nn8l : lEver et al.ll2ni5l ). Ground-based RVs with <1 km s“^ precisions will there¬ 


fore continue to play a critical role in confirming this nearby population of young stars in the 
foreseeable future. 
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Reference s. — B13=lBo wler et alj ll2013lV B15 = [ B owler et al.l ||201^: D14=|Di ttmann et al.l 1120141 ): CIO = | Chauvin e t al.| |2010[): GOlj^ray e^a^ l|2^ ^); H 96= [Hawlev et al 


'^9€); J12=|jansoi^^al.l |2(^ll): Lll1 201ll ): L13jLe£ine_^^ 20131: MlS^Malo et al.||201 j; Mel4 ^Me R04= Reid et al 


' 20od) : R07= lReid et alJ 1I2OO7I); R08 ^Roser et al.l ll2008l) : RlO^Roeser et al.| ll2010l) : S98= IStauffer et alJ lll9^ : Sh09= IShkolnik et al.l ll2009l) : Shl2^ 
TW=This work: Z05 jzacharias et al.l 1I2OO5I) : Z13=UCAC4 llZacharia^^il.l 1I2OI3I) ) 


200^|_^;0^ ^a^^^aL, 

Shkol^k et al.l ||2012^ : 
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Table 2. Companion Properties 


Name 

Ang. Sep 
(") 

Proj. Sep 
(AU) 

SpT 

YMG? 

Age 

(Myr) 

Mass 

(ATj.p) 

Discovery 

Reference 

2MASS J01225093-2439505 B 

1.5 

52 ± 6 

L5 ± 1 

AB Dor? 

120 ± 10 

13-25 

B13 

2MASS J02155892-0929121 C 

3.5 

90 ± 10 

M7 ± 1 

Tuc-Hor? 

30-40 

42 ± 15 

BIO, TW 

2MASS J02594789-0913109 B 

0.6 

22 ± 8 



>4000 

87 ± 13 

J12, TW 

IRXS J034231.8+121622 B 

0.8 

19.1 ± 0.9 

LO ± 1 


60-300 

35 ± 8 

B15 

HD 23514 B 

2.6 

360 ± 3 

M8 ± 1 

Pleiades 

120 ± 10 

60 ± 10 

R12 

G 160-54 C 

3.3 

70 ± 30 

M7.0 ± 0.5 


>1000 

85 ± 16 

B15 

2MASS J05464932-0757427 B 

2.8 

130 



>1000 

87 ± 3 

J12 

2MASS J06475229-2523304 B 

1.1 

260 ± 40 

M3.0 ± 0.5 




TW 

2MASS J08540240-3051366 B 

1.7 

17 ± 7 

M5.0 ± 0.5 


<200? 

~40 ± 20 

CIO 

2MASS J11240434+3808108 B 

8.3 

170 ± 11 

M9.5 ±0.5 

UMa? 

500 ± 100 

81 ± 5 

C03 

PYC J11519+0731 B 

0.5 

18 ± 3 

M8 ± 1 


>10 

15-90 

TW 

G 180-11 B 

1.6 

45 ± 5 

M6.0 ± 0.5 


<200 

90-150 

MOl, D07 

2MASS J15594729+4403595 B 

5.6 

150 ± 11 

M7.5 ± 0.5 


50-200 

43 ± 9 

J12, B15 

GJ 4378 Ab 

0.5 

7 ± 1 

M8 ± 1 


>4000 

90 ± 3 

TW 


^Estimated spectral type based on colors. 


References. — B10= |Bergfo rs et al.l ll2O10| b B13=|^ wiere t al.l ll2 0lj); B15=[ j owler e t al.l ll2Q15 [l: CQ3=|C lose et ^ 
1 I 2 OO 3 ); C10^ Chau^n et^.| 12010^ : D07 jpaemgen et al.l J200^ J12= jjanson et ^ <|2012|) : MOl ^McCarthv et al.l il20^K : 
R12= Rodriguez et al .1 ll201i li: TW=This work 






























Table 3. Adaptive Optics Imaging Observations 


Name 

Date 

(UT) 

Epoch 

(UT) 

Telescope/ 

Instrument 

Filter 

N X Coadds x texp 

(s) 

Separation 

(mas) 

P.A. 

n 

A mag 

(mag) 

m.B*’ 

(mag) 

2M01225093-2439505 AB 

2013 Aug 17 

2013.6258 

Keck II/NIRC2 

Y 

14 X 20 X 2 

1448 ± 3 

216.52 ± 0.09 

7.6 

± 

0.3 

10.61 ± 0.03 

18.2 ± 0.3 

2M02155892-0929121 AB 

2012 Aug 23 

2012.6451 

Keck II/NIRC2 

Ks 

20 X 100 X 0.028 

575.5 ± 0.5 

289.98 ± 0.03 

1.17 

± 

0.02 

7.86 ± 0.02 

9.03 ± 0.03 


2012 Oct 12 

2012.7816 

Subaru/IRCS 

K 

5 X 5 X 0.5 

575 ± 2 

289.9 ± 0.5 

1.16 

± 

0.02 

7.84 ± 0.03'^ 

9.00 ± 0.03® 


2013 Aug 17 

2013.6259 

Keck II/NIRC2 

Ks 

10 X 10 X 0.2 

576.35 ± 0.11 

289.351 ± 0.011 

1.112 

± 

0.006 

7.88 ± 0.02 

8.99 ± 0.02 

2M02155892-0929121 AC 

2012 Aug 23 

2012.6451 

Keck II/NIRC2 

R^+corbOO 

5 X 1 X 2.9 

3430 ± 3 

299.4 ± 0.2 

4.1 

± 

0.2 

7.57 ± 0.02 

11.7 ± 0.2 


2012 Oct 12 

2012.7816 

Subaru/IRCS 

K 

4 X 5 X 0.5 

3370 ± 80 

299.57 ± 0.17 

4.50 

± 

0.16 

7.54 ± 0.03'^ 

12.04 ± 0.16® 

2M02155892-0929121 Aab 

2012 Aug 23 

2012.6451 

Keck II/NIRC2 

Ks 

20 X 100 X 0.028 

42 ± 7 

112 ± 2 

1.2 

± 

0.3 

7.86 ± 0.08 

9.1 ± 0.2 


2013 Aug 17 

2013.6259 

Keck II/NIRC2 

Ks 

10 X 10 X 0.2 

42 ± 7 

308 ± 5 

1.17 

± 

0.05 

7.86 ± 0.03 

9.03 ± 0.04 

2M02594789-0913109 AB 

2013 Aug 17 

2013.6259 

Keck II/NIRC2 

Ks 

10 X 10 X 0.2 

648.6 ± 0.3 

320.53 ± 0.01 

2.70 

± 

0.02 

10.22 ± 0.02 

12.92 ± 0.03 

IRXS J034231.8+121622 AB 

2007 Dec 14 

2007.9513 

Keck II/NIRC2 

Ks 

10 X 30 X 1.5 

883.0 ± 0.2 

17.58 ± 0.09 

3.62 

± 

0.04 

9.313 ± 0.018 

12.93 ± 0.04 


2010 Aug 29 

2010.6585 

Gemini-S/NICI 

H 

8 X 1 X 60 

851 ± 3 

18.7 ± 0.1 

3.84 

± 

0.06 

9.580 ± 0.019 

13.42 ± 0.06 

2M05464932-0757427 AB 

2013 Jan 17 

2013.0448 

Keck II/NIRC2 

Ks 

10 X 10 X 1.0 

2877.3 ± 1.2 

114.48 ± 0.03 

3.86 

± 

0.03 

9.81 ± 0.02 

13.67 ± 0.04 

2M06475229-2523304 AB 

2012 Oct 13 

2012.7847 

Subaru/IRCS 

J 

5 X 5 X 0.4 

1116 ± 8 

28.9 ± 0.5 

5.22 

± 

0.09 

8.36 ± 0.02 

13.58 ± 0.09 


2012 Oct 13 

2012.7847 

Subaru/IRCS 

H 

4 X 1 X 10.0 

1071 ± 7 

28.7 ± 0.3 

4.89 

± 

0.12 

7.77 ± 0.04 1 

12.66 ± 0.13 


2012 Oct 13 

2012.7847 

Subaru/IRCS 

K 

4 X 10 X 0.26 

1089 ± 23 

28.8 ± 0.2 

5.1 

± 

0.5 

7.53 ± 0.02S:::. 

12.6 ± 0.5® 


2013 Jan 18 

2013.0478 

Keck II/NIRC2 

K 

5 X 100 X 0.106 

1086.6 ± 1.2 

28.08 ± 0.07 

5.538 

± 

0.012 

7.53 ± 0.02^ 

13.07 ± 0.03® 


2014 Nov 8 

2014.8531 

Keck II/NIRC2 

Ks 

9 X 10 X 0.1 

1088 ± 2 

28.10 ± 0.09 

5.33 

± 

0.16 

7.55 ± 0.02 

12.88 ± 0.16 


2014 Nov 8 

2014.8531 

Keck II/NIRC2 

i^s+cor600 

10 X 1 X 30 

1089.9 ± 0.7 

27.80 ± 0.03 






2M08540240-3051366 AB 

2013 Jan 18 

2013.0478 

Keck II/NIRC2 

j 

4 X 100 X 0.053 

1741 ± 3 

159.08 ± 0.03 

3.00 

± 

0.07 

9.08 ± 0.03 

12.08 ± 0.07 


2013 Jan 18 

2013.0478 

Keck II/NIRC2 

H 

4 X 100 X 0.053 

1741.2 ± 1.3 

159.07 ± 0.08 

2.97 

± 

0.03 

8.46 ± 0.04 

11.43 ± 0.05 


2013 Jan 18 

2013.0478 

Keck II/NIRC2 

K 

4 X 100 X 0.181 

1740.9 ± 1.8 

159.07 ± 0.05 

2.848 

± 

0.008 

8.142 ± 0.03'^ 

10.99 ± 0.03® 

2M11240434+3808108 AB 

2014 Jun 13 

2014.4463 

P60/Robo-AO 

i' 

1 X 1 X 60 

8250 ± 90 

131.2 ± 0.5 

4.63 

± 

0.05 

12.554 ± 0.001^ 

17.18 ± 0.05^ 


2014 Jun 13 

2014.4463 

P60/Robo-AO 

z' 

1 X 1 X 60 

8270 ± 90 

131.1 ± 0.5 

3.90 

± 

0.03 

13.985 ± 0.013^ 

17.88 ± 0.03^ 

PYC 11519+0731 AabB 

2012 May 22 

2012.3896 

Keck II/NIRC2 

Y 

11 X 5 X 0.5 

496 ± 3 

110.6 ± 0.3 

5.3 

± 

0.3 

9.31 ± 0.03‘i 

14.6 ± 0.3*^ 


2012 May 22 

2012.3896 

Keck II/NIRC2 

J 

20 X 100 X 0.05 

497 ± 2 

110.4 ± 0.2 

5.40 

± 

0.15 

8.82 ± 0.03 

14.22 ± 0.15 


2012 May 22 

2012.3896 

Keck II/NIRC2 

H 

22 X 10 X 0.028 

496.7 ± 1.4 

110.4 ± 0.2 

5.4 

± 

0.2 

8.14 ± 0.04 

13.5 ± 0.2 


2012 May 22 

2012.3896 

Keck II/NIRC2 

Ks 

17 X 100 X 0.028 

497.1 ± 1.1 

110.46 ± 0.09 

5.23 

± 

0.09 

7.90 ± 0.03 

13.13 ± 0.10 


2012 May 22 

2012.3896 

Keck II/NIRC2 

V 

10 X 10 X 0.11 

500 ± 3 

110.4 ± 0.2 

4.97 

± 

0.13 

7.80 ± 0.11® 

12.77 ± 0.17® 


2012 Jun 24 

2012.4798 

Keck II/NIRC2 

H 

20 X 100 X 0.028 

498.2 ± 2.3 

110.21 ± 0.15 

5.18 

± 

0.09 

8.14 ± 0.04 

13.32 ± 0.10 


2012 Jun 24 

2012.4798 

Keck II/NIRC2 

Ks 

19 X 100 X 0.028 

499.2 ± 0.5 

110.30 ± 0.10 

5.08 

± 

0.12 

7.90 ± 0.05 

12.98 ± 0.13 


2013 Jan 18 

2013.0479 

Keck II/NIRC2 

i7+cor600 

20 X 1 X 30 

500.9 ± 1.8 

108.76 ± 0.05 







2013 Apr 19 

2013.2957 

Gemini-S/NICI 

H 

10 X 10 X 1.14 

517 ± 10 

108.3 ± 0.5 

6.0 

± 

0.2 

8.14 ± 0.04 

14.1 ± 0.2 


2013 Apr 19 

2013.2957 

Gemini-S/NICI 

Ks 

10 X 10 X 1.14 

511 ± 10 

107.9 ± 0.5 

5.4 

± 

0.2 

7.89 ± 0.04 

13.3 ± 0.2 


2013 Dec 18 

2013.9628 

Keck II/NIRC2 

Rs+corbOO 

9 X 2 X 10 

514 ± 2 

107.15 ± 0.04 











Table 3—Continued 


Name 

Date 

(UT) 


Epoch 

(UT) 

Telescope/ 

Instrument 

Filter 

N X Coadds x £exp 

(s) 

Separation 

(mas) 

P.A. 

(°) 

Amag 

(mag) 

(mag) 

G 180-11 AB 

2013 Feb 

04 

2013.0949 

Keck II/NIRC2 

Y 

10 X 10 X 0.5 

1620.4 ± 1.3 

254.98 ± 0.03 

2.089 

± 

0.008 

9.60 ± 0.02'^ 

11.70 ± 

0.02'^ 


2013 Feb 

04 

2013.0949 

Keck II/NIRC2 

J 

10 X 10 X 0.2 

1619.9 ± 1.2 

254.92 ± 0.05 

2.063 

± 

0.010 

9.08 ± 0.02 

11.14 ± 

0.02 


2013 Feb 

04 

2013.0949 

Keck II/NIRC2 

H 

19 X 10 X 0.028 

1619.6 ± 1.4 

254.97 ± 0.03 

2.08 

± 

0.02 

8.41 ± 0.03 

10.49 ± 

0.04 


2013 Feb 

04 

2013.0949 

Keck II/NIRC2 

Ks 

15 X 10 X 0.028 

1620.5 ± 1.3 

254.93 ± 0.04 

2.06 

± 

0.03 

8.19 ± 0.03 

10.25 ± 

0.04 


2014 Jun 

13 

2014.4463 

P60/Robo-AO 

r' 

1 X 1 X 60 

1610 

± 

50 

253.2 ± 1.6 

2.18 

± 

0.03 

13.544 ± 0.004^ 

15.724 ± 

0.027f 


2014 Jun 

13 

2014.4463 

P60/Robo-AO 

i' 

1 X 1 X 60 

1640 

± 

50 

253.9 ± 1.6 

2.09 

± 

0.02 

11.749 ± 0.003f 

13.839 ± 

0.017*' 


2014 Jun 

13 

2014.4463 

P60/Robo-AO 

z' 

1 X 1 X 60 

1630 

± 

50 

253.8 ± 1.6 

2.03 

± 

0.02 

11.142 ± 0.005^ 

13.172 ± 

0.018* 

2M15594729+4403595 AB 

2014 Jun 

13 

2014.4463 

P60/Robo-AO 

i' 

1 X 1 X 60 

5670 

± 

70 

284.4 ± 0.6 

6.37 

± 

0.07 

10.433 ± 0.001^ 

16.80 ± 

0.07* 


2014 Jun 

13 

2014.4463 

P60/Robo-AO 

z' 

1 X 1 X 60 





5.55 

± 

0.05 

10.147 ± O.OOlf 

15.70 ± 

0.05* 

GJ 4378 Aab 

2012 Oct 

13 

2012.7842 

Subaru/IRCS 

J 

4 X 1 X 30.0 

473.4 

± 

1.2 

348.20 ± 0.18 

2.907 

± 

0.011 

8.71 ± 0.02 

11.62 ± 

0.02 


2012 Oct 

13 

2012.7842 

Subaru/IRCS 

H 

5 X 1 X 8.0 

475.6 

± 

1.8 

348.3 ± 0.16 

3.14 

± 

0.07 

8.13 ± 0.03 

11.27 ± 

0.07 


2012 Oct 

13 

2012.7842 

Subaru/IRCS 

K 

5 X 30 X 0.145 

475.3 

± 

1.7 

348.4 ± 0.10 

3.02 

± 

0.03 

7.93 ± 0.05'= 

10.95 ± 

0.06" 


2013 Jan 

18 

2013.0471 

Keck II/NIRC2 

K 

13 X 10 X 0.053 

492.8 

± 

1.5 

347.85 ± 0.13 

3.06 

± 

0.03 

7.92 ± 0.05" 

10.98 ± 

0.06" 


2013 Aug 

17 

2013.6258 

Keck II/NIRC2 

Y 

10 X 5 X 0.5 

526.7 

± 

0.4 

347.22 ± 0.03 

3.353 

± 

0.014 

9.21 ± 0.02'^ 

12.56|± 

0.03'* 


2013 Aug 

17 

2013.6258 

Keck II/NIRC2 

J 

10 X 10 X 0.05 

526.0 

± 

0.3 

347.22 ± 0.03 

3.11 

± 

0.03 

8.70 ± 0.02 

ll.Slj^ 

0.04 


2013 Aug 

17 

2013.6258 

Keck II/NIRC2 

H 

10 X 10 X 0.05 

526.1 

± 

0.2 

347.18 ± 0.03 

3.107 

± 

0.019 

8.13 ± 0.03 

11.24^ 

0.03 


2013 Aug 

17 

2013.6258 

Keck II/NIRC2 

Ks 

15 X 10 X 0.1 

525.96 

± 

0.16 

347.15 ± 0.03 

3.066 

± 

0.008 

7.87 ± 0.03 

10.94'± 

0.03 


2013 Aug 

17 

2013.6258 

Keck II/NIRC2 

K 

10 X 10 X 0.1 

526.06 

± 

0.08 

347.194 ± 0.018 

3.046 

± 

0.018 

7.93 ± 0.05" 

10.97 ± 

0.05" 


2013 Aug 

17 

2013.6258 

Keck II/NIRC2 

V 

20 X 10 X 0.114 

525.9 

± 

0.6 

347.23 ± 0.04 

2.819 

± 

0.011 

7.58 ± 0.10*= 

10.40 ± 

0.11" 


2014 Aug 

04 

2014.5901 

Keck II/NIRC2 

Ks 

20 X 10 X 0.2 

577.2 

± 

0.4 

346.2 ± 0.02 

3.06 

± 

0.03 

7.87 ± 0.03 

10.93 ± 

0.04 


^Magnitude of the primary (rriA) decomposed from the integrated-light magnitude of the A and B components (mA n) as follows: mA = ^ab + 2.5 log(l -|- 10 ^™a.g/2.5^ 
Uncertainties are derived in a Monte Carlo fashion. Unless otherwise noted, J, and Ks magnitudes are from 2MASS JSkrutskie et al.]l2Q06ll . 

'^Magnitude of the companion (ms) is computed from Amag and mA- 

‘^Integrated-light Xg-band magnitude is converted to Kmko using relations from [Leggett et~ 1 1200^ . 

■^y-band integrated-light magnitude is derived from the typical Y — J color of dwarfs for the system spectral type from iRavner et al.l ll2009ll and the integrated-light J-band 
magnitude. 

^L^band integrated-light magnitude is derived from the typical Ks — L' color of dwarfs for the system spectral type from Ic^limowski et al. I ll2004l~) and the integrated-light 
X^-band magnitude. 

^Unresolved i'^ and z' magnitudes are from the Sloan Digital Sky Survey DR9 lAhn et al.ll20lj) . Note that optical variability of the primary and/or companion can affect the 
decomposed photometry of one or both components. 




















Table 4. Spectroscopic Observations 


Object 

Date 

(UT) 

Telescope/ 

Instrument 

Filter 

Slit Width 

(") 

Plate Scale 
(mas pix~^) 

Exp. Time 
(min) 

Resolution 

(=AA/A) 

Standard® 

2MASS J01225093-2439505 A 

2014 Jun 28 

SOAR/Goodman 

GG455 

0.46 


5 

1800 

LTT 6248 

2MASS J02155892-0929121 A 

2012 Dec 28 

Keck I/HIRES 


0.86 


1 

45000 

GJ 908 


2014 Jun 28 

SOAR/Goodman 

GG455 

0.46 


4 

1800 

LTT 6248 


2014 Jun 28 

SOAR/Goodman 

GG455 

0.46 


5 

5900 

LTT 6248 


2014 Dec 08 

CTIO 1.5-m/CHIRON 




45 

28000 

HD 42581 

2MASS J02155892-0929121 C 

2014 Jan 18 

IRTF/SpeX-SXD 


0.5 

150 

12 

1200 

HD 13936 

2MASS J02594789-0913109 A 

2012 Dec 28 

Keck I/HIRES 


0.86 


5 

45000 

GJ 908 

HD 23514 B 

2013 Feb 3 

Keck I/OSIRIS 

Hbb 


50 

45 

3800 

HD 23258 


2013 Feb 3 

Keck I/OSIRIS 

Kbb 


50 

42 

3800 

HD 23258 

G 160-54 A 

2007 Dec 19 

CFHT/ESPaDOnS 


0.53 


16.7 

68000 

GJ 273 


2012 Jan 29 

UH 2.2-m/SNIFS 



400 

7 

1300 

HR 1544, G191B2B 

G 160-54 G 

2012 Sep 24 

IRTF/SpeX-SXD 


0.5 

150 

24 

1200 

HD 25792 

2MASS J06475229-2523304 A 

2013 Apr 30 

du Pont/Echelle 


0.75 


10 

45000 

GJ 273, GJ 433 


2014 Dec 3 

UH 2.2-m/SNIFS 



400 

3.3 

1300 

Feige 67 


2014 Dec 10 

CTIO 1.5-m/CHIRON 




15 

28000 

HD 217357 

2MASS J06475229-2523304 B 

2014 Dec 8 

Keck I/OSIRIS 

Kbb 


20 

18 

3800 

HD 58886 

2MASS J08540240-3051366 A 

2009 Jun 13 

du Pont/Echelle 


0.75 


15 

45000 

GJ 908 

2MASS J08540240-3051366 B 

2014 Dec 8 

Keck I/OSIRIS 

Hbb 


20 

10 

3800 

HD 58886 


2014 Dec 8 

Keck I/OSIRIS 

Kbb 


20 

16 

3800 

HD 58886 

2MASS J11240434+3808108 B 

2011 Jan 22 

IRTF/SpeX-SXD 


0.5 

150 

8 

1200 

HD 105388 

PYC J11519+0731 Aab 

2012 Jun 1 

IRTF/SpeX-SXD 


0.8 

150 

3 

750 

HD 97585 


2012 Jul 5 

Keck I/HIRES 


0.9 


3.6 

4500 



2013 Feb 22 

ESO-MPG 2.2-m/FEROS 




20 

48000 



2014 Jun 9 

CFHT/ESPaDOnS 


0.53 


6.7 

68000 

GJ 625 


2014 May 22 

Mayall/RC-Spectrograph 

GG495 

1.5 

690 

4 

2500 

HZ 44 

PYC J11519+0731 B 

2012 Jun 25 

Keck I/OSIRIS 

Kbb 


20 

40 

3800 

HD 116960 

G 180-11 B 

2013 Feb 1 

Keck I/OSIRIS 

Jbb 


50 

6 

3800 

q Her 


2013 Feb 1 

Keck I/OSIRIS 

Hbb 


50 

6 

3800 

q Her 


2013 Feb 1 

Keck I/OSIRIS 

Kbb 


50 

6 

3800 

q Her 

2MASS J15594729+4403595 A 

2012 Feb 12 

CFHT/ESPaDOnS 


0.53 


5 

68000 

GJ 628 


2012 Aug 23 

UH 2.2-m/SNIFS 



400 

1.7 

1300 

HR5501, EG131 

2MASS J15594729+4403595 B 

2012 Sept 9 

Keck II/ESI 


1.0 

150 

40 

4000 

G191B2B, Feige 67 

GJ 4378 A 

2009 Aug 20 

du Pont/Echelle 


0.75 


10 

45000 

GJ 908, GJ 699 


2014 Jun 28 

SOAR/Goodman 

GG455 

0.46 


4 

1800 

LTT 6248 












Table 4—Continued 


Object 

Date 

(UT) 

Telescope/ 

Instrument 

Filter 

Slit Width 

(") 

Plate Scale 
(mas pix~^) 

Exp. Time 
(min) 

Resolution 

(=AA/A) 

Standard^ 


2014 Jun 28 

SOAR/Goodman 

GG455 

0.46 


5 

5900 

LTT 6248 

GJ 4378 Ab 

2013 Jul 31 

Keck I/OSIRIS 

Hbb 


50 

28.5 

3800 

HD 219833 


2013 Jul 31 

Keck I/OSIRIS 

Kbb 


50 

11 

3800 

HD 219833 

GJ 4379 B 

2009 Aug 20 

du Pont/Echelle 


0.75 


10 

45000 

GJ 908, GJ 699 


2014 Jun 28 

SOAR/Goodman 

GG455 

0.46 


4 

1800 

LTT 6248 


2014 Jun 28 

SOAR/Goodman 

GG455 

0.46 


5 

5900 

LTT 6248 


^For our high-resolution optical spectra, this refers to the RV standard. For our near-infrared and low-resolution optical 
spectra, this refers to telluric standards. 
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Table 5. Ha Emission and Optical Spectral Classification 


Name Ho EW^ Ti05 CaHl CaH2 CaH3 Adopted SpT 



(A) 

Index 

Index 

Index 

Index 

(±0.5) 

2MASS J01225093-2439505 A 

-5.0 

0.417 

0.783 

0.424 

0.677 

M4.0 

2MASS J02155892-0929121 Aab 

-5.6 

0.482 

0.810 

0.485 

0.716 

M3.5 

G 160-54 Aab 

-0.8 

0.413 

0.820 

0.416 

0.701 

M4.5 

2MASS J06475229-2523304 Aab 

-0.3 

0.890 

1.009 

0.950 

0.952 

K7.0 

PYC J11519+0731 Aab 

-2.8 

0.594 

0.866 

0.612 

0.814 

M2.0 

2MASS J15594729+4403595 A 

-2.3 

0.685 

0.877 

0.648 

0.821 

M1.5 

2MASS J15594729+4403595 B 

-6.4 

0.195 

0.828 

0.303 

0.674 

M7.5 

GJ 4378 Aal2 

-4.5 

0.393 

0.769 

0.408 

0.670 

M4.5 

GJ 4379 B 

0.1 

0.457 

0.801 

0.459 

0.725 

M3.5 


^Negative values indicate emission. 
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Table 6. Radial Velocities 


Name 

Component 

UT Date 
(YYYY-MM-DD) 

RV 

km s“^ 

Reference 

2MASS J01225093-2439505 

A 

2010-09-24 

11.4 ± 0.2 

M14 


A 

2011-08-24 

12.3 ± 2.5 

K13 


A 

2012-12-28 

9.6 ± 0.7 

B13 

2MASS J02155892-0929121 

Aa 

2007-11-08 

0.6 ± 5.0 

K13 


Aa 

2010-07-26 

0.5 ± 0.3 

M14 


Aa 

2010-09-19 

-0.6 ± 0.3 

M14 


Aa 

2012-02-11 

8.3 ± 0.3 

M14 


Aa 

2012-07-18 

10.1 ± 0.6 

K14 


Aa 

2012-12-28 

6.2 ± 0.4 

TW 


Aa 

2014-12-08 

8.3 ± 0.6 

TW 

2MASS J02594789-0913109 

A 

2012-12-28 

4.9 ± 0.5 

TW 

IRXS J034231.8+121622 

A 

2005-12-21 

35.4 ± 0.4 

S12 

2MASS J06475229-2523304 

Aa 


-56.5 

T06 


Aa 

2009-02-28 

-67 ± 6 

K13 


Aa 

2011-09-14 

-24.4 ± 0.8 

M14 


Aa 

2013-04-30 

-2.1 ± 0.9 

TW 


Aa 

2014-12-10 

-5.5 ± 0.8 

TW 

2MASS J08540240-3051366 

A 

2009-06-13 

44.5 ± 0.6 

TW 

2MASS J11240434+3808108 

A 

2006-05-11 

-11.5 ± 0.5 

S12 


B 


-14 ± 3 

R09 

PYC J11519+0731 

Aa 

2012-07-05 

50.9 ± 1.0 

TW 


Aa 

2013-02-22 

14.7 ± 1.0 

TW 


Aa 

2014-06-09 

69.5 ± 1.0 

TW 


Ab 

2012-07-05 

-51.9 ± 1.0 

TW 


Ab 

2013-02-22 

-18.2 ± 1.0 

TW 


Ab 

2014-06-09 

-75.5 ± 1.0 

TW 


7 


-0.3 ± 1.0 

TW 

G 180-11 

A 

2006-08-13 

-15.5 ± 0.7 

S12 

2MASS J15594729+4403595 

A 

2013-03-03 

-15.8 ± 0.5 

M14 


A 

2012-02-12 

-19.6 ± 0.6 

TW 

GJ 4378 A 

al 

2009-08-20 

-7.0 ± 0.8 

TW 


a2 

2009-08-20 

49.2 ± 1.1 

TW 

GJ 4379 B 


2009-08-20 

20.6 ± 0.4 

TW 


References. 


Male et alJ 1 

20141); 

Torres et al. 

120061) 


TW = this work. 
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Table 7. IRXS J034231.8+121622 AB Astrometry 


Epoch 

Separation 

(mas) 

PA 

(°) 

Ref 

2007.95 

883.0 ± 0.2 

17.58 ± 0.09 

TW 

2008.63 

860 ± 8 

17.3 ± 0.4 

J12 

2008.87 

866 ± 8 

17.8 ± 0.4 

J12 

2010.66 

851 ± 3 

18.7 ± 0.1 

TW 

2012.02 

834 ± 57 

17.6 ± 1.7 

J14 

2012.65 

831 ± 2 

18.71 ± 0.07 

B15 

2013.04 

822 ± 8 

19.1 ± 0.7 

B15 


Ref erences. — B 15=|Bowler et all 
1201^: J12^Janson et all 1 I 2 OI 2 I) : 

J14= |janson et al.l ||2014| ): TW=this work. 
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Table 8. Allers & Liu 2013 Gravity Indices 


Parameter 

2MASS J02155892-0929121 C 

G 160-54 G 

2MASS J11240434-I-3808108 B 

VO 2 Index 

1.023 ± 0.011 

1.042 ± 0.007 

1.069 ± 0.011 

FeH 2 Index 

1.082 ± 0.012 

1.109 ± 0.010 

1.210 ± 0.015 

K Ij Index 

1.030 ± 0.006 

1.072 ± 0.005 

1.111 ± 0.007 

if-cont Index 

1.005 ± 0.006 

0.968 ± 0.005 

0.931 ± 0.007 

FeHj Index 

1.066 ± 0.010 

1.113 ± 0.010 

1.202 ± 0.019 

EW(Na I 1.138 /rm) (A) 

6.9 ± 0.7 

12.1 ± 0.2 


EW{K I 1.169 iim) (A) 

2.7 ± 0.5 

4.3 ± 0.4 

5.5 ± 0.5 

EW(K I 1.177 iim) (A) 

4.0 ± 0.5 

6.6 ± 0.4 

8.8 ± 0.5 

EwIk I 1.253 ^J,m) (A) 

2.4 ± 0.4 

4.5 ± 0.4 

5.4 ± 0.4 

Gravity Score^ 

lnl2 

OnOO 

OnOO 

Gravity Class 

INT-G 

FLD-G 

FLD-G 


^Gravity scores are for the FeH, VO, alkali lines, and ff-band continuum shape. See I Allers &: Tii] ll2013l~) 
for details. 
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Fig. 1.— Low-resolution SOAR/Goodman (i?~1800) and U H2.2ni/SNIFS (R^1300 ) optical spec¬ 
tra of host stars to ultracool companions. Red templates from iBochanski et al.l (1200711 show the bes t 
matches based on visual and index-based classification with Hammer (Table O ICovev et al.l 120071 1. 
Half-integer subtypes are created by averaging two templates together. All spectra are normalized 
between 7200-7400 A and offset by a constant. Althongh these spectra originate from instruments 
with different resolving powers and wavelength spans, our spectral types agree with those in the 
literature (to within ±0.5 subtypes) for nearly all stars which have been previously classified. 
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Fig. 2.— Keck/NIRC2 R-band image of 2MASS J01225093-2439505 B. Contours represent 0.03%, 
0.5%, 1%, 10%, and 50% of the peak flux after convolution with a Gaussian kernel with a FWHM 
equal to that of the image PSF. North is up and east is left. 
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X(k) 


Fig. 3.— Moderate-resolution (i?~5900) spectra of 2MASS J02155892-0929121 A, GJ 4378 A, and 
GJ 4379 B from the SOAR/Goodman spectrograph. None of the targets show signs of Li I A6708 
absorption. All spectra are normalized between 7200-7400 A and offset by a constant. 
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Fig. 4.— Keck/NIRC2 iiT^-band coadded image of the young 2MASS J02155892~0929121 system. 
The primary (“Aab”) is positioned behind the 600 mas partly-opaque coronagraph. The stellar 
(“B”) and brown dwarf (“C”) companions are located at separations Cf.T and 3'.'5, respectively. 
North is up and east is left. 
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Epoch Epoch 



Fig. 5.— Astrometry of the substellar companion 2MASS J02155892"0929121 C between 2008- 
2014 relative to the primary, which itself is a close 40 mas visual binary. The solid curve shows 
the expected relative astrometry of a background object; the gray shaded regions show 1 cr and 
2 cr confidence intervals. Open symbols represent the predicted astromet ry for a stationary ob - 
ject at the epochs of the observations. Compared to astrometry from iBergfors et ahl (j2010l l. 
2MASS 402155892-0929121 C is clearly comoving. There is some evidence of orbital motion in 
the system, which is likely caused by motion from the Aab primary, Aab-B component, and/or 
between the Aab-C components of this quadruple system. 
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Fig. 6.— NIRC2 images of 2MASS J02155892-0929121 A and B in 2012 and 2013. The “A” 
component is slightly elongated compared to “B” and is well-fit as a marginally-resolved visual 
binary utilizing the PSF of “B”. The bottom two rows from left to right show the raw image of 
“A,” the residuals after jointly fitting a two-star PSF of “B” and subtracting the the brighter 
component (“Aa”), the same after subtracting the fainter component (“Ab”), and the residuals 
after subtracting both. In August 2012, 2MASS J02155892-0929121 Ab is ~40 mas southeast of 
Aa, and in August 2013 it is ~40 mas northwest of Aa. 
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Fig. 7.— IRTF/SpeX 0.8-2.4 jum spectrum of 2MA SS J02155892-0929121 C (black). Visual 
and index-based classification from lAllers &: Liul (|20l3l 1 yields a spectral type of M7 ± 1 and an 
intermediate surface gravity (“INT-G”). A reduced surface gravity is evident from the shallow alkali 
lines and FeH band streng th in 2MASS J02155892-0929121 C compared to field templates (red) 
from ICushing et al.l (j2005l i and iRavner et al.l ([20091) . This implies the companion is quite young 
(30-150 Myr) and probably substellar, in line with the system’s likely association with the Tuc-Hor, 
/3 Pic, or Columba moving groups. All templates are Gaussian smoothed to match the resolving 
power of our SpeX spectrum (R~1200). All spectra are normalized and offset by a constant. 
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Fig. 8.— Orbit fits for the intermediate-age (60-300 Myr) M4-I-L0 binary IRXS J034231.8-1-121622 
AB. Although the astrometric baseline (2007-2013) spans a small fraction of the orbital period, 
only a subset of Keplerian orbits are consistent with the observations. Each curve represents a draw 
from the posterior distribution of allowed orbital elements, conditioned on the astrometric data. 
The orbits are plotted only from the expected position of the companion in 1980 to the expected 
position in 2050. In the upper panel, the black dot denotes the position of the primary star. The 
lower panel shows a magnihed version of the upper panel. 
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Fig. 9.— Posterior probability distributions for the eccentricity and total mass of 

IRXS J034231.8+121622 AB (see Figure [8]). Initial orbit fits of this system point to a low ec¬ 
centricity (<0.3) and a total mass of RiO.1-0.2 Mq. 
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Fig. 10.— Keck/OSIRIS 1.4-2.4 /xm spectrum of the substellar companion to the dusty F6 Pleiades 
star HD 23514. HD 23514 B best resembles field M7-L0 templates (red) from the I RTF Spectral Li¬ 


brary. Our OSIRIS H and K spectra (black) are flux calibrated to photometry from lRodrieuez et al 


(120121 1. All spectra are Gaussian smoothed to a common resolving power, normalized, and offset 
by a constant. 
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Fig. 11.— IRTF /SpeX 0.8-2.4 /xm spectrum of G 160-54 C (black). Visual and index-based 
classification from lAllers &: Liul (j2013l i yields a near-infrared spectral type of M7 ± 0.5 and a field 
surface gravity (FLD-G). Field templates (red) are from the IRTF Spectral Library and have been 
Gaussian smoothed to match the resolving power of our SpeX spectrum (i2~1200). All spectra are 
normalized and offset by a constant. 



























Fig. 12.— UH 2.2-m/SNIFS spectru m of 2MASS J0647 5229 -2523304 A (black) com pared to field 
K- and M-dwarf templates (red) from iMann et al.l (|2013l l and iBochanski et al.l (|2007l l. respectively. 
2MASS J06475229-2523304 A broadly resembles the K7V/M0V templates but the detailed molec¬ 
ular band strengths are not a good match. Instead, index measurements point to a low surface 
gravity. The Ha line is slightly in emission {EW=-0.3 A). Spectral regions with strong telluric 
absorption have been removed and all spectra are normalized and offset by a constant. 
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Fig. 13.— 4"x4" Subaru/IRCS and Keck/NIRC2 AO images of 2MASS J06475229-2523304 AB 
between 2012 and 2013. The companion is easily resolved, with a separation of T.'l at both epochs. 
Contours represent 0.02%, 0.5%, 1%, 5%, and 50% of the peak flux after convolution with a 
Gaussian kernel with a FWHM equal to that of the image PSF. North is up and east is left. 
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Epoch 


Fig. 14.— Relative astrometry of 2MASS J06475229-2523304 AB between 2012 and 2014 (filled 
symbols). The solid curve shows the expected path of a background object and gray shaded regions 
show 1 a and 2 a confidence intervals. Open symbols represent the predicted astrometry for a 
stationary object at the epochs of the observations. The pair are clearly comoving and physically 
bound. 
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Fig. 15.— Keck/OSIRIS K-band spectrum of 2MASS J06475229-2523304 B (black). Despite being 
5.5 mag fainter than the K7 host star 2MASS J06475229-2523304 A in the NIR, this companio n 

i2nn9l i. 


has a spectral type of M3V based on comparison to field template (red) from iRavner et al. 

This suggests one or both components are not on the main sequence; the most likely scenario is the 
system is a distant (240 pc) K7 subgiant with an M3 compaion. Our OSIRIS spectrum has been 
Gaussian smoothed to match the resolving power of the templates (R~2000), and all spectra are 
normalized and offset by a constant. 










~ 70 - 


. ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' . 

■ HgO FeH H 2 O 

:— 1 —'—'—'— 1 —'—'—'— 1 —: 

. HgO H 2 O . 


: M3V : 

M4V : 




- IV16V ^ ” 


■ M7V ^ ■ 



: • i?=2ooo : 

: •i?= 2 ooo : 



Fig. 16.— Keck/OSIRIS 1.5-2.4 spectrum of 2MASS J08540240-3051366 B (black). The 
best match is to M4-M6 field templates (red) from the IRTF Spectral Library; we adopt a near- 
infrared spectral type of M5 ± 0.5. Our OSIRIS spectrum has been Gaussian smoothed to match 
the resolving power of the templates (R~2000), and all spectra are normalized and offset by a 
constant. 
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Fig. 17.— IRTF/SpeX 0.8-2.4 um spectrum of 2M ASS J11240434+3808108 B (black). Visual and 
index-based classification from lAllers &: Liul (j2013l l yields a near-infrared spectral type of M9.5 ± 0.5 
and a field snrface gravity (FLD-G). Field templates (red) are from the IRTF Spectral Library and 
the M6V-M9V objects have been Gaussian smoothed to match the resolving power of our SpeX 
spectrum (i?~1200). The LO template has i2~750 so our science spectrum has been smoothed to 
match that resolution. All spectra are normalized and offset by a constant. 
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Fig. 18.— Optical P60/Robo-AO images of G 180-11 AB, 2MASS J11240434-I-3808108 AB, and 
2MASS J15594729-1-4403595 AB. All three late-M companions are easily detected in 60-s frames. 
North is up and east is left. 
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Fig. 19.— Mayall/RC-Spectrograph medium-resolution (i2Ri2500) optical spectrum of 

PYC J11519-I-0731 A (black). Visual classification yields a spect ral type of M2, t houg h some 
peculiarities are evident compared to SDSS templates (red) from iBochanski et al.l (|2007l i. The 


inset shows the Ha emission line, which is marginally resolved here but well-separated in our high- 
resolution spectrum. A joint ht to both components with the AMOEBA downhill-simplex algorithm 
yields a differential velocity of ~154 km s“^. However, for our RV analysis in this work we only 
include the high-resolution spectra. All spectra are normalized and offset by a constant. 
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Fig. 20.— IRTF/SpeX SXD unresolved spectrum of PyC J1 1519+0731 Aab (black) compared to 
the M2.0V template G1 806 (red) from iRavner et al.l (|2009l b Several absorption features in our 
spectrum are unusually deep such as the Mg I Al.503/1.505 ^m and Si I A1.589 /rm lines, possibly 
reflecting a higher metallicity. The template is smoothed to the resolution of our spectrum (R~750) 
and each individual band is normalized. 
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Fig. 21.— Top: The rest position of the Li I A6708 line in PYC J11519+0731 Aab (solid gray) 
and predicted positions of the same line in the primary and secondary given their measured RVs 
(dashed red) at three epochs. We do not find convincing evidence for lithium absorption in either 
component. Bottom: Similar to the top panels but for Ha emission lines in the primary and 
secondary components. The emission lines are easily split in our 2012 and 2014 data but have a 
small relative RV in our 2013 data set. 
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Fig. 22.— Radial velocities of PYC J11519+0731 Aa and Ab. Our three epochs enable a mea¬ 
surement of the mass ratio (q) and systemic velocity ( 7 ) for the system. The negative slope of the 
best-fit line (solid line) is equal to the mass ratio of the system (0.95 ± 0.01), and we measure a 
systemic RV of-0.3 ± 1.0 km s“^. 
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Fig. 23.— Partially constrained UVW kinematics and XYZ space positi o ns o f PYC 
J11519+0731 Aab compared to young moving group members from iTorres et al.l (120081 ). Red 
circles show the stars’ position for distances of 20, 30, 40, 50, and 60 pc based on the systemic RV. 
At close distances of Ri20-30 pc, PYC J11519+0731 Aab is consistent with {3 Pic moving group, 
but physically would be tens of parsecs from known members. 
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Fig. 24.— l"5xl"5 Keck/NIRC2 R-, J-, H-, Ks-, and L'-band AO images of PYC J11519+0731. 
The (f!5 companion PYC J11519+073 B is clearly visible in these unsaturated coadded frames. 
Contours represent 0.3%, 0.6%, 1%, 10%, and 50% of the peak flux after convolution with a 
Gaussian kernel with a FWHM equal to that of the image PSF. North is up and east is left. 
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Fig. 25.— NIRC2 i7-band coronagraphic observations of PYC J11519+0731 AabB. The Aab com¬ 
ponent, unresolved here, is centered behind the ~0.1%-throughput 600 mas diameter coronagraph. 
The image on the left is a single frame in our ADI sequence; PYC J11519+0731 B (A77=5.2 mag) 
is visible but resembles speckle noise. The right panel shows the result after PSF subtraction with 
a library of NIRC2 PSF reference images. North is up and east is left. 
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Fig. 26.— Relative astrometry of PYC J11519+0731 Aab and B between 2012 and 2014 (filled 
symbols). The solid curve shows the expected path of a background object and gray shaded 
regions show 1 a and 2 a confidence intervals. Open symbols represent the predicted astrometry 
for a stationary object at the epochs of the observations. The pair is clearly comoving and exhibits 
slight but significant orbital motion. 
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Fig. 27.— Keck/OSIRIS iT-band spectrum of PYC J11519+0731 B (black). Compared to field 
templates (red) from the IRTF Spectral Library (left panel), PYC J11519+0731 B best resembles 
the M8 object so we adopt a near-infrared spectral ty pe of M8 Y 1. Th e right panel shows a 
comparison to “very low-gravity” templates (blue) from lAllers fc LujI ( 20131 ). none of which match 
our spectrum as well as the high-gravity counterparts. Our spectrum has been Gaussian smoothed 
to match the resolving power of the comparison object. All spectra are normalized and offset by a 
constant. 
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Fig. 28.— 7 a contrast curve (left) and companion sensitivity map (right) for the PYC J11519+0731 
system. Our limiting contrast of AT/sslS mag at 1" corresponds to planetary masses (<13 Mjup) at 
young ages (<1 Gyr) and masses in the brown dwarf regime for ages beyond ~1 Gyr. The conversion 


from contrast to planet mass assumes hot s tart ( Cond) evolutionary models from iBaraffe et al 
([2003) and circular orbits. See iBowler et al.l (|2015l l for details. 
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Fig. 29.— Keck/OSIRIS 1.2-2.4 //m spectrum of G 180-11 B (black). The best match is to M5-M6 
held templates (red) from the IRTF Spectral Library. We adopt a near-infrared spectral type of 
M6 ± 0.5. Our OSIRIS spectrum has been Gaussian smoothed to match the resolving power of 
the templates (R~2000) and all spectra are normalized and offset by a constant. 



























- 84 - 


I ' ' ' ' I ' ' ' ' I ' ' ' ' I 
2MASS J15594729+4403595 B 


Ha 

I Tip 

Tio caH rn 
r I-1 


■n 


TiO 


r 

©Op VO Na I 


M6.5V[ 



6000 6500 7000 7500 8000 8500 

X(A) 

Fig. 30.— Keck/ESI optical spectrum of 2MASS J15594729+440 3595 B (black). The best matches 
are M7.0 and M7.5 templates (red) from iBochanski et al.l (|2007l h where the half-types are formed 
from the average of neighboring integer templates. We adopt an optical spectral type of M7.5 ± 0.5 
but note that the ~7400 A VO feature appears weaker than the field template. 
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Fig. 31.— Complete 0.5-2.4 /im spectrum of 2MASS J15594729+4403595 B (black) compared to 
the high-gravity M8 field object LP 412-31 (red) from the IRTF Spectral Library. The optical 
(0.4-0.9 /rm) region is our new Keck/ESI spectrum while the near-infrared (0.9-2.4 ^m) region is 
our IRTF/SpeX spectrum from iBowler et al.l (|2015l L The lithium absorption and angular iL-band 
shape in 2MASS J15594729-I-4403595 B point to an age of Ri50-200 Myr, making this a nearby 
young system that does not appear to belong to any known moving groups. 
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Fig. 32.— Relative astrometry of GJ 4378 Aa and Ab between 2012 and 2015. The solid curve shows 
the expected path of a background object and gray shaded regions show 1 a and 2 a confidence 
intervals. The pair are clearly comoving and exhibit significant orbital motion. 
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Fig. 33.— l75xl75 Keck/NIRC2 Y-^ J-, H-, Ks--, and L'-band AO images of GJ 4378 Aab. 
Contours represent 0.3%, 0.6%, 1%, 10%, and 50% of the peak flux after convolution with a 
Gaussian kernel with a FWHM equal to that of the image PSF. North is up and east is left. 
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Fig. 34.— Keck/OSIRIS 1.5-2.4 ^um spectrum of GJ 4378 Ab (black). The best match is to M8-M9 
field templates (red) from the IRTF Spectral Library. We adopt a near-infrared spectral type of 
M8 ± 1. Note that the slightly angular iL-band shape, which points to a young age of <150 Myr, is 
at odds with the lack of Ha emission in the wide comoving companion GJ 4379 B, which suggests an 
old age of several Gyr. Our OSIRIS spectrum has been Gaussian smoothed to match the resolving 
power of the templates (i?~750-2000) and all spectra are normalized and offset by a constant. 
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Fig. 35.— The known ultracool companions to stars. Systems analyzed in this work are denoted 
with red squares and mostly have late-M spectral types. The dearth of young (<100 Myr) com¬ 
panions later than ~L5 reflects the paucity of young planets disc o vered via direct im aging. Note 


that we have excluded the HR 8799 pla nets (Marois et al. 20081 : iMarois et al.N2010l l. GJ 504 b 


(jKuzuhara et al.ll2013l l. and HD 95086 b (jRameau et al.ll2013l l because their spectral types are ei¬ 
ther poorly cons trained or may defy conventional classification schemes. Known companions (grey 
circles) are from iDeacon et al.l (120141 1 and are supplemented with our own compilation. 












































































